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We present a detailed analysis of classical solutions in the bosonic sector of the 
electroweak theory which describe vortices carrying a constant electric current X. 
These vortices exist for any value of the Higgs boson mass and for any weak mix- 
ing angle, and in the zero current limit they reduce to Z strings. Their current 
is produced by the condensate of vector W bosons and typically it can attain bil- 
lions of Amperes. For large 2 the vortices show a compact condensate core of size 
~ embedded into a region of size ~ I where the electroweak gauge symmetry is 
completely restored, followed by a transition zone where the Higgs field interpolates 
between the symmetric and broken phases. Outside this zone the fields are the same 
as for the ordinary electric wire. An asymptotic approximation of the large I solu- 
tions suggests that the current can be arbitrarily large, due to the scale invariance 
of the vector boson condensate. Finite vortex segments whose length grows with I 
seem to be perturbatively stable. This suggests that they can transfer electric charge 
between different regions of space, similarly to thunderbolts. It is also possible that 
they can form loops stabilized by the centrifugal force - electroweak vortons. 
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I. INTRODUCTION 

More than 20 years ago Witten proposed a field theory model that admits classical 
solutions describing vortices carrying a constant current - 'superconducting strings' [l|. 
This model has a local U(l)xU(l) invariance and consists of two copies of the Abelian 
Higgs model for fields {A^i^\(j)i) and (y4|f'',02)- It also includes an interaction between the 
two complex scalars chosen such that in vacuum one has 0i 7^ but 02 = so that the vector 
field A^^^ is massive while A^f^ is massless and can be identified with the electromagnetic 
field. 

The model admits as a solution the Abrikosov-Nielsen-Olesen (ANO) vortex made 
of the 'vortex fields' (yl|(^\0i), with vanishing 'condensate fields' (y4|f'',02)- This embedded 
vortex is however unstable, but being topological it does not unwind into vacuum and relaxes 
to a 'dressed vortex' which contains a condensate of charged scalar bosons in the core and 
has 02 7^ 0. Giving then a non-trivial phase to the condensate field 02 produces a current 
and promotes the 'dressed' vortex to the superconducting string supporting the long-range 
Biot-Savart field represented by A^f^ 7^ 0. 



The Witten string superconductivity has been much studied j^, Q], mainly in the cos- 
mological context [sl, j^, since Witten's model can be viewed as sector of some high energy 
Grand Unification Theory (GUT) jl| that could perhaps be relevant at the early stages of 
the cosmological evolutions. Using the typical values of the GUT parameters for estimates 
gives for the string current enormous values of order 10^° Amperes, which suggests interest- 
■„g applicafous g Str.ng superconduct.vity ■„ the GUT-related „o,.AbeMa„ models 
has also been studied [7[, in which case the string current is produced by a condensate of 
charged vector bosons. 

Although the GUT physics could be important, one may wonder whether a similar string 
superconductivity could exist also in a less exotic context, at lower energies, as for example 
in the electroweak sector of Standard Model. In fact, the U(l)xU(l) symmetry of Wit- 
ten's model is contained in the SU(2)xU(l) electroweak gauge symmetry. In addition, the 
electroweak theory contains a pair of complex Higgs scalars, one of which could well be 
responsible for the formation of the vortex while the other one - for the condensate. Since 
;he ANO vortices can be embedded into the electroweak theory in the form of Z strings 
Sj, one might expect that current-carrying generalizations of the latter could exist. These 
would be electroweak analogues of Witten's superconducting strings. 

However, no attempts to construct such solutions have ever been undertaken. This can 
probably be explained by the following reasons. The current-carrying Witten strings are 
usually viewed as excitations over the 'dressed' currentless vortex obtained by minimizing 
the energy of the 'bare' embedded ANO vortex. This explains their essential properties, as 
for example the value of the critical current Sj. Now, the 'bare' electroweak Z strings are 



also unstable [9|, and it was conjectured 
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ll| that they could similarly relax to 'W-dressed 



Z strings'. However, a systematic search for such solutions gave no result p^]. This can 
probably be explained by the fact that Z strings are non-topological and can unwind into 
vacuum [l3^, so that there is little chance that they could be stabilized by the condensate. 
But if one does not find currentless 'dressed' Z strings, it is natural to think that current- 
carrying strings in the electroweak theory do not exist either. This presumably explains why 
there has been almost no activity on electroweak vortices during the last 10-15 years. 

At the same time, it is difficult to believe that Z strings are the only possible vortex 
solutions in the electroweak theory. The theory admits rather non-trivial solutions, such 
as sphalerons 14\, periodic BPS solutions [isl, spinning dumbbells [l^, oscillons jisl ]. 



spinning sphalerons [19|], and others (see [20| for a review). In addition, in the semi 
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hmit where the SU(2) field decouples the theory admits current-carrying vortices [211], so 
that it is plausible that they could exist also for generic values of the weak mixing angle. 

We have therefore analysed the problem and found that superconducting vortices indeed 
exist in the electroweak theory, despite the non-existence of the 'W-dressed Z strings'. In 
other words, these two types of solutions are not necessarily related. To construct the 
solutions, we essentially reverse the standard 'engineering' procedure used in Witten's model. 
There one starts from the 'dressed' currentless vortex and increases it 'winding number 
density' that determines the phase of the condensate; in what follows we shall call this 
parameter 'twist'. This produces a current that first increases with the twist but then starts 
to quench and finally vanishes when the solution reduces to the 'bare' ANO vortex The 
superconducting strings thus comprise a one-parameter family that interpolates between the 
'dressed' vortex and the 'bare' ANO vortex. 

We construct this family in the opposite direction, by starting from the 'bare' vortex and 
then decreasing its twist. Within Witten's model this gives of course the same solutions but 
in the reversed order - their current first increases then starts to quench and vanishes for 
zero twist when the solution reduces to the 'dressed' string. However, the advantage of our 
method is that it can be applied also within the Weinberg-Salam theory, where there are no 
'dressed' currentless vortices but only the 'bare' ones - Z strings. Their twist is determined 
by the eigenvalue of the second variation of the energy functional. Decreasing the twist 
gives us solutions with a non-zero current. Further decreasing it shows that the current 
always grows and tends to infinity when the twist approaches zero. As a result, contrary to 
what one would normally expect, and presumably because of the vector character of their 
current-carriers, the electroweak vortices do not (generically) exhibit the current quenching, 
so that they do not need to admit the 'dressed' currentless limit. 

The stability analysis of these vortices shows that their finite segments are perturbatively 
stable. More precisely, so far this property has been demonstrated only in the semilocal 
limit 221], but it is very likely that the result extends for generic values of the weak mixing 
angle. The length of the stable segments increases with the current and can in principle 
attain any value, since there is no upper bound for the current. 

This may have interesting consequences. First of all, this suggests that loops made of 
the stable segments and balanced by the centrifugal force could perhaps be stable as well. 
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Although studying this goes far beyond the scope of the present paper, the possibihty to 
have stable solitonic objects in the Standard Model could be very important, since if such 
electroweak vortons exist, they could perhaps contribute to the dark matter. 

Another, perhaps a more direct consequence, is related to the well-know fact that, since 
Z strings are non-topological, they can have finite length. This suggests that their current- 
carrying generalizations could also exist in the form of finite (and stable) segments connecting 
oppositely charged regions of space. They would then be similar to thunderbolts. For the 
solutions we could explicitly construct the current can be as large as 10^° Amperes, which 
exceeds by several orders of magnitudes the power of the strongest thunderbolts in the 
Earth atmosphere. This suggests that superconducting vortices could be important for the 
dynamics of the Standard Model, although analysing this issue in detail would again lead 
us too far away from our main subject. 

The present paper is devoted to a systematic analysis of the string/ vortex- type solutions 
in the electroweak theory (terms string and vortex are assumed to be synonyms). In what 
follows we shall show that every Z string admits a three-parameter family of non-Abelian, 
current carrying generalizations. In the gauge where the radial components of the gauge 
fields vanish, the upper and lower components of the Higgs field doublet behave quite anal- 
ogously to the vortex field and condensate field in Witten's model. The new vortices exist 
for (almost) any value of the weak mixing angle and for any Higgs boson mass. They show 
a compact, regular core containing the W-condensate that carries a constant electric cur- 
rent, producing the Biot-Savart electromagnetic field outside the vortex. In the comoving 
reference frame the vortex is characterized by the current I and by the electromagnetic 
and Z fluxes through its cross section. After performing a Lorentz boost, it develops also 
a non-zero electric charge density, as well as momentum and angular momentum directed 
along the vortex. 

It seems that the vortex current can be as large as one wants, at least we could not find an 
upper bound for it. We have also managed to construct a simple approximation in the large 
current limit which suggests that the current can be arbitrarily large, due to the fact that 
it is carried by vectors whose condensate exhibits the scale invariance. To the best of our 
knowledge, a similar effect has never been reported before. For large currents the charged 
W boson condensate is contained in the vortex core of size ~ 1/X, which is surrounded by 
a large region of size ~ I where the Higgs field is driven to zero by the string magnetic 
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field so that the electroweak gauge symmetry is completely restored. However, this does not 
destroy the vector boson superconductivity, since the scalar Higgs field is not the relevant 
order parameter in this case. Outside the symmetric phase region there is a transition zone 
where the Higgs field relaxes to the broken phase and the theory reduces to the ordinary 
electromagnetism. This reminds somewhat of the electroweak vacuum polarization scenario 
discussed by Ambjorn and Olesen |15t], in which the system stays in the Higgs vacuum if 
only the magnetic field is not too strong, while a very strong field restores the full gauge 
symmetry. 

The rest of the paper is made maximally self-contained. In Sec. II the essential elements of 
the Weinberg-Salam theory are introduced. Sec. HI presents the reduction to the stationary 
sector, the derivation of the conserved quantities, and the discussion of the Bogomol'nyi- 
Prasad-Sommerfield (BPS) limit. Sec. IV contains the further reduction to the axially sym- 
metric sector, the basic field ansatz (14.41) . and the field equations (I4.12p - (l4.20p . In Sec.V the 
boundary conditions are discussed, while Sec. VI and Sec. VII contain the derivation of the 
charges and a brief description of the known solutions. The new solutions are first considered 
in Sec. VIII in the small current limit, when they can be treated as small deformations of Z 
strings. In Sec. IX they are presented at the full non perturbative level for generic values of 
the current. Sec.X describes the large current limit. Solutions for special parameter values 
are considered in Sec. XI, while Sec. XII contains concluding remarks. 

There are also three appendices. Appendices A and B contain the derivation of the local 
solutions at the symmetry axis and at infinity. Appendix C describes the superconducting 
strings in Witten's model. This Appendix is made self-consistent, but it uses the same 
methods and notation as in the main text, so that it illustrates our procedure on a simple 
example. 

A very brief and preliminary summary of our results has been announced in Ref. 23 1. 

Recently there has been an intense activity on the non-Abelian strings with gauge group 
SU(N)xU(l) 24]. These strings have nothing to do with ours, since they are obtained 
within the context of a supersymmetric theory that is different from the electroweak theory 
even for N=2. They are currentless and are supposed to be relevant in the QCD context, 
as models of gluon tubes. There has also been a report on current-carrying strings 25|], but 
again in a model with a different Higgs sector in which solutions taking values in the Cartan 
subalgebra of SU(N)xU(l) are possible. 
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II. WEINBERG-SALAM THEORY 



The bosonic part of the Weinberg-Salam (WS) theory is described by the action 
Here the complex Higgs field $ is in the fundamental representation of SU(2) and 



(2.1) 



B^v = dfjBj, — d^Bfj, , 



=\^d,-^B,-^ T^W;j (2.2) 

where t"' are the Pauli matrices and c in (12.11) is the speed of light. The SU(2) gauge field 
and the U(l) field B^, are sometimes called isospin and hypercharge fields, respectively. 
Here and below we denote dimensionfuU quantities by boldfaced symbols. In order to pass 
to dimensionless variables, let us first introduce the dimensionless coupling constants 

5' = — = cos6'w, 5'' = — = sin^w, (2.3) 
go go 

where the physical value of the Weinberg angle is 

^'2 = sin2^w = 0.23... (2.4) 

The parameter 



go = VFTF (2.5) 

is related to the electron charge via e = gg'hcgo so that its dimension is [go] = [1/e] and 
the value 

{gg'gor = (2.6) 



AttHc An ° ' 137 
The Higgs field vacuum expectation value $o has the dimension [$o] = [^/L] and the value 

$0 = 54.26 X 10^ Volts. (2.7) 

Introducing the length scale L = l/go$o = 1-52 x 10~^^ metres (-\/2L is the Z boson 
Compton length) one can define the dimensionless variables 

B W ^ I (t>i\ 8A 

x'^ = go$ox^ B, = g'^, W;: = ^,-^, <^' = ^=r , /3 = -T, (2-8) 

*o *o ^0 \ 6r, So 



9 



such that the action (12. ip becomes 



S = ^ / Cd^x. (2.9) 



ego 



Here 



and 



= "i^ W^.W'^^^ - ^ fi^.fi^'^ + (D^$)tD^$ - ^ ($t$ - 1)' , (2.10) 



D,<t>=(^d,-'-B,-'-T'^W^^^. (2.11) 
The Lagrangian (I2.10p is invariant under the SU(2)xU(l) gauge transformations 

$ ^ U$, W ^ UWU"^ + 2iVdf,V-^dx^ , (2.12) 

with 

2 2 



U = exp ( Je + ^r"r ) (2.13) 



where and 6'" are functions of x'^ and 

W = {B^ + T''W''^)dx'' (2.14) 

is the SU(2)xU(l) Lie-algebra valued gauge field. Varying the action with respect to the 
fields gives the field equations, 

d^'B,, = ((D.<l>)t$ - $tD.$) ^ (2.15) 
D'-y^U = 9^\ ((D,<I>)V'^$ - <|.tr'^D,$) ^ <^V,^ (2.16) 
+ ^ ($^$ - 1)$ = 0, (2.17) 

with DfjW'^p = d^W^p + eabc^^^ap- Varying the action with respect to the spacetime 
metric gives the energy momentum tensor, 

dC 

T''u = '^9'"'g^-K^, (2.18) 

which evaluates to 

T"^ = -l^W'^'^'Wl^ - -^B'^'B^^ + (D^$)tD^$ + (D^$)"^D'^$ - S^^C . (2.19) 
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Let us define vacuum as 

/ 1 \ 

(2.20) 

Any field with T^^ = is a gauge transformation of tlie vacuum. Let us consider small 
fluctuations around the vacuum in the unitary gauge, assuming that WJ^ and -B^ are small 

and that $ = with real and small. Linearising the field equations fl2.15p -( l2.17p 

V " / 

with respect to the perturbations gives 

d^F^" = 0, 
d^Z'"' + 77112" = 0, 
5MWr + ^wW: = (a = 1,2), 

d^,d>'^ + = 0. (2.21) 



Here the electromagnetic and Z fields are defined as 



= - ^ Z,, = B,, + Wj„ (2.22) 



with the potentials 



A^ = j^B^-^Wl, Z^ = B, + Wl, (2.23) 



while the field masses 



mz = mw = gm^, = v^mz. (2.24) 

v2 

Multiplying these by hcg^^Q = e^o/{gg') gives the dimensionfuU masses, for example 

mzc^ = e^o/iV2gg'), (2.25) 



which evaluates 91.18 GeV [26|]. The value of the parameter /3 defining the Higgs field 
mass is currently unknown, however, there are indications that it belongs to the interval 
1.5 < /? < 3.5 [261. 



The electromagnetic and Z fields in (12.221) are defined under the assumption that the 
system is close to the vacuum. In what follows we shall need to define these fields also off 
the vacuum. This can be done in a gauge invariant way, however there is no unique definition 
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in this case [27j], although all known definitions reduce to fl2.22p when the fields approach 



the vacuum. Unless the otherwise is stated, we shall adopt the definition of Nambu 



16| 



F^. = ^B^,-^ n-W;, , Z,, = B,, + n^YJ% , (2.26) 

where 

n^ = ^^T''^, ^ = <^/V¥^. (2.27) 

It seems that this definition gives the most sensible results as compared to the other defini- 
tions. Its only disadvantage is that the 2-forms (12.261) are not closed in general, so that there 
are no field potentials. However, this cannot be considered as a drawback, since there is 
no reason why the Maxwell equations should hold off the Higgs vacuum. It is worth noting 
that a similar preference in favour of the Nambu definition was made also in Ref. {29 ]. 



Another well known definition of the fields is according to 't Hooft 
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where 



with 



-^/^l - ^ - ^Vfiu, Z^u = + Vfiu (2.28) 



V^. = n'^W;, - eabcn^D^n'D.n' = d^V^ - d,V^ (2.29) 



= n'^Wl + 2iS}d^i , = di^n^ + Ca^cW^^ (2.30) 

The field strengths are closed in this case and the potentials are known explicitly. In the 
unitary gauge where n" = ^3 the potentials reduce simply to A,^ = — ^ and 

Z^ = B^ + in which form they are most often used (see, e.g., [15:]). However, as we shall 
see below, the t' Hooft definition gives sometimes singular results, while applying instead 
the Nambu definition makes things perfectly sensible. 

III. SYMMETRY REDUCTION 

In what follows we shall be considering stationary and translationally symmetric fields 

invariant under the action of two spacetime isometrics generated by the Killing vectors 

d d 
Kio) = Q^, ^(3) = (3.1) 

Since these vectors commute and all the internal symmetries of the theory are gauged, there 



exists a gauge where the symmetric fields do not depend on x° = t and = z 



3Q| 
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Let (Jq, be a constant (co)vector in the plane and 5"q, = ea/sc^ its dual. These two 

vectors are orthogonal and have the norm 

= -a'^a^ = = -(^o)' + (a^)^ . (3.2) 

Using these two vectors, the fields respecting the symmetries (13.11) can be written as 



>V =[uaa + uaa]dx'^ + Bidx' + r"[M„(T„ + + r^W^dx', $ = , (3.3) 




where a = 0,3 and i,k = 1,2, while m, m, -u", -Bj, W°, 0i, 02 depend on x , so that there 
are in total 20 independent real functions. The ansatz ( 13. 3p keeps its form under gauge 
transformations (12.121) with O, 9"" depending only on x^. 

Inserting this ansatz into (12.101) gives the reduced Lagrangian 



+ ^ 



- \Dk^f - ^ $t(^ + UaT^y^ - ^^(n + Mar")2<l.^ - ^ (<Dt$ - 1)2 (3.4) 

where DkUa = dkUa + eabc^luc and = (9^ — | Bk — | r"W^)$. Varying this gives the 
field equations, where, as is evident from the structure of the reduced Lagrangian, one can 
consistently set 

u = u'' = 0. (3.5) 
The remaining non-trivial equations then read 

dAu = ^^^iu + Uar'')<^, (3.6a) 

D,DsUa = ^^^iUa + UT'')<^, (3.6b) 

dsBsk = I ^i^^Dk^ - (Dfc$)t$) , (3.6c) 

DsW:, - ahabcUtDkU, = I ^(<l>tr-D,<l> - (D^-K) V»<1>) , (3.6d) 

DsDs $- — (« + UaT^f^ = ^ ($^$ - 1)$ . (3.6e) 

Let us discuss the global quantities associated to solutions of these equations - the magnetic 
fiuxes, charge and current, energy and momentum. 
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A. Fluxes 

The electromagnetic and Z fields produce magnetic fluxes through the plane, 

*. = i/e..f„<i'x. = (3.7) 

The dimensionfull values are = ^ f e^'^Fjfcrf^x = "^p/gQ and = ^z/go- The fluxes 
are gauge invariant, but their values depend on the definition of F^^^ and Z^^, chosen. If 
the fields are defined according to 't Hooft, then they admit potentials, in which case one 
can use the Stokes theorem to express the fluxes in terms of integrals over the boundary at 
infinity. With the Nambu definition one should calculate the surface integrals. 

In general there is no reason for the fluxes to be quantized. Indeed, let us briefly recall 
the standard argument in favour of the vortex flux quantization. Far away from the vortex 
the condition 



{d.-'-B,-'- r»W;) $ = (3.8) 



should be fulfilled, whose solution is 



$(x) = Pexp (^^ j{B^ + T''YJl)dx^^ <l>(xo) . (3.9) 

Here V stands for the path ordering along a contour C interpolating between xq and x. 
Choosing C to be a large circle around the vortex gives the holonomy condition 

$(xo) =Pexp (^-j{B^ + T''W;)dx^^^{xo). (3.10) 

If = 0, as for example in the semilocal limit (7 — > (see Sec JXI Al below), then the 
condition (I3.10p reduces to 

]- ! B^dip = 27cn (3.11) 
2 Jo 

and so the hypercharge field flux is quantized, \E'_b = 47m. 

Let us now suppose that = = but 7^ and let us denote = J^^ W<^ dip. 
The condition (13.101) then gives 

01 (xo) = et(*-+*-)0i(xo), M^o) = e^(*--*-V2(xo) , (3.12) 



and if, for example, 4'2{xo) = 0, then one should have "^b + "^w = 47rn but the difference 
"^B — "^w can be arbitrary. The fluxes are therefore not quantized in general. 
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B. Current 

The conserved electromagnetic current is 

J, = 9^F^,. (3.13) 

Since nothing depends on = (x°, x^), one has Ja = —OkF^a- This being a total derivative, 
integrating over the plane gives 

Ia = J Jad^x = - j dkFka dPx = - j) CksFka dx'' . (3.14) 

Here /q is the total electric charge per unit length of the vortex and is the total electric 
current through the x^x"^ plane. Since the integration in (13.141) is performed over the bound- 
ary at infinity, in the region where all fields approach the vacuum and the electromagnetic 
field is uniquely defined, the value of la does not depend on the chosen definition of 
Let us consider the dimensionfull version of (13.131) . 

^J, = 5^F^„ (3.15) 

where F^j, = go^'g-^^i' ^^^^ the derivative is taken with respect to = a;'^/go$o- One has 

j J„rf^X=-C j dkFkad^X= -C^o j dkFka d^ X = C^ola, (3.16) 

from where it follows that la gives the current in units of 

C$0 = c X 54.26 X 10^ Vohs = 1.8 x 10^ Amperes. (3.17) 

This value is quite large. Below we shall present vortex solutions whose current is typically 
J3 ~ 1 — 10, which looks modest but corresponds in fact to billions of Amperes ! Very 
large currents are typical for superconducting strings. In the GUT-related model of Witten 
the current can be as large as 10^° Amperes - because the GUT Higgs field vacuum 

expectation value is much larger than $o- 

C. Energy and momentum 

Contracting the Killing vectors with the energy-momentum tensor (I2.19P gives conserved 
Noether currents, 

J^j,^=T^,K'^ 5,j5,) = 0. (3.18) 
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Integrating over the plane gives the conserved charges per unit vortex length - the 

vortex energy and momentum densities, 

E = [ TlKr^.d^x = [ T%d^x, (3.19) 



P = / TlKY^.d^x = / T%d^x . (3.20) 



The dimensionfuU values are E = $q and P = $g P/c. 

Let us consider the momentum. According to Ref. 32|, in gauge field theories without 
global internal symmetries the density j^j^^ for a spacelike Killing vectors K has the total 
derivative structure. Calculating the energy-momentum tensor ( 12.19P for the field ( 13. 3p 
shows that T\ is indeed a total derivative, 

9 9 

= CToCTs SfcSfc (3.21) 

with 

9 9 2 

= ^^'^^ + ^ . (3-22) 

where the field equations fl3.6ap - fl3.6el) have been used to pass to the second line. For globally 
regular fields the total momentum is therefore completely determined by the boundary term, 

P = ao(T3 y efcsSfcrfx* . (3.23) 

Let us now consider the energy. It is given by a sum of two manifestly positive terms 
which we shall call, respectively, electric and magnetic energy, 

E = Jt^ d^x = £2) d^x =Ei + E2. (3.24) 

The structure of Si is very similar to that of T3 in (13.211) . 

2 I 2 

£1 = -^^—^ — - dk^ki (3.25) 
so that it is a total derivative. The magnetic energy is 

^2 = ^ (W^,)^ + -^,{B,kf + \DM' + ^m'- (3.26) 
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D. Bogomol'nyi bound and the electroweak condensate 

It is instructive to mention the existence of the non-trivial lower bound for the magnetic 
energy, observed by Ambjorn and Olesen [isl. Using the identity 

= i I A$ + ^e^,D,^^ + i ^^$^(5,^ + r'^W^fc)* - ieskds{<^^ D^^) (3.27) 



the density £2 can be rewritten as 

£2 = ■^Awl + ^t^k^^T''^ 



2 J Ag'^ 



2 



,2 



+ ^ eife5,fc - ie,kdsi<^^Dk<^) - ^ . (3.28) 

From this it follows that for (3 > 1 one has 

E2> Ji^eikBik-teskdsi^^Dk^)-^)cfx, (3.29) 

which defines the lower energy bound. This bound is, however, somewhat special. Consid- 
ering just one vortex and integrating over the whole plane, the first term in the integrand 
in fl3.29p gives the magnetic flux, the second term gives zero, while the third one gives the 
infinite area of the plane with a negative coefficient. The conclusion is that E2 is always 
larger than minus infinity, which is true but trivial. 

A less trivial conclusion can be made if one considers an infinite number of parallel 
vortices forming a periodic lattice in the plane, since the field Bik then does not 

vanish at infinity and its contribution can overcome the effect of the area term. It is then 
sufficient to integrate just over the elementary periodicity cell with the area A. One uses the 
relation Bik = gg' Fik + g'^Zik where Zik does not contribute to the flux since its potentials 
(defined in according to t' Hooft) are gauge invariant and therefore periodic if Z^k is 



periodic [l5|. The bound then reads 

where is the flux of Fik through A. If the system contains only a constant magnetic field 
F12 > m^/e then this relation gives 

E2 = \ {F,2fA > ^(Fi2 - ^)A > 0. (3.31) 
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This shows that for large F12 the system becomes unstable, since its energy grows as (-^12)^ 
while the lower energy bound grows only as F12. The bound can be achieved for (3 = 1 ii 
the first order Bogomol'nyi equations are fulfilled, 

Wrfc + yQfc$tr'^<l> = 0, (3.32a) 
+ =0, (3.32b) 

A$ + ieijDj<^ = . (3.32c) 

Periodic solutions of these equations describe a non-linear condensate of W,Z and Higgs 
fields forming an infinite number of non-Abelian vortices that spontaneously appear in the 



very strong magnetic field [15|. It turns out that non-trivial solutions exist if only 



15| 



.2 ^2 



— < {F12) < — (3.33) 
e e 

where (-F12) is the magnetic field averaged over A. If the lower inequality here is violated, 
then the Higgs field falls to the vacuum, while in the opposite limit it is driven to zero fisl . 
This can be interpreted as the electroweak symmetry restoration in the very strong magnetic 
field. 

Since the magnetic fields satisfying the conditions f l3.33p are extremely strong, it was 



suggested to use the Witten superconducting strings of the GUT origin as their source |15 |. 
At certain distances from the string core the magnetic field falls into the interval (13.331) . so 
that an electroweak condensate can form within a cylindrical shell around the string. Such 
systems are sometimes called 'W-dressed superconducting strings' [l^, {31 ]. 

Below we shall construct electroweak vortices showing a compact superconducting core 
embedded into a symmetric phase region that is wrapped in a cylindrical shell where the 
Higgs field interpolates between the symmetric and broken phases. This reminds very much 
of the 'W-dressed superconducting strings', up to the fact that our strings are of purely 
electroweak origin and have nothing to do with the GUT physics. They can exist for any 
value of f3, which is why the Bogomol'nyi equations are of little use for us. 

The above arguments go differently in the semilocal limit, where g = (see Sec JXI All . 
since the area term in (I3.30p then vanishes and the magnetic flux through the whole plane 
is no longer constrained to be infinite. The flux is then quantized as = 47rn and the 
energy lower bound is 27m. Eq. (l3.32ap implies then that W";, = so that the Yang-Mills 
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field can be gauged away, while the remai ning Bogomol'nyi equations (13.32bp . fl3.32cp admit 



vortex solutions for any given G Z 



35], |36| 



IV. AXIAL SYMMETRY 

Let us now further restrict our consideration to fields which are invariant with respect 
to rotations in the plane. This means that in addition to the two Killing vectors 

(13. ip there exists a third one. Choosing the polar coordinates in the plane, = pcosip, 
x"^ = p sin yj, it is given by 



^(^) = — • (4-1) 



d_ 

Since this vector commutes with the other two, for fields respecting all three symmetries 
at the same time there exists a gauge where they do not depend neither on nor on (p. 
They are then given by Eg. (13.31) with all the functions depending only on p. The field 
equations (13.61) become in this case a system of 16 ordinary differential equations for the 16 
independent functions left in the ansatz (13.31) after imposing the condition (13. 5p . 

Within the gauge chosen one can still perform gauge transformations (I2.12p with B, 6°- 
depending only on p. We use this freedom to set to zero the radial components of the gauge 
fields, 

Bp = 0, W^ = 0. (4.2) 

This reduces the number of independent functions to 12. Now, the inspection of the field 
equations reveals that one can consistently set to zero the imaginary components of the 
fields. These include the projection of the Yang-Mills field and also the imaginary part 
of the Higgs field. We can therefore set on-shell 

Wj = 0, $5(<1>) = 0. (4.3) 

As a result, we are left with only 8 independent real functions, in which case we can 
parametrize the remaining non-trivial components in (13. 3p as 

yy = u{p) ffadx" — v{p) dip + [ui{p) ffadx" — vi{p) dip] 

+ [uz{p) a^dx'^ - vsip) dip], $ = f ^'^^^ , (4.4) 

V Up) J 

where a = 0, 3 and /i, /2 are real. This ansatz has the following properties. 
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a. It is invariant under spacetime symmetries generated by 

d d d 



b. It is invariant under complex conjugation, 



(4.5) 



W = W*, $ = $*. (4.6) 

c. It keeps its form under Lorentz rotations in the x^ plane, whose effect on (14. 4p is 
to Lorentz-transform the components of the (co) vector cTq: 

ctq o"o cosh a — (73 sinh a, cosh a — ctq sinh a, (4-7) 

where a is the boost parameter. The norm cr^ = cr| — ctq is Lorentz-invariant, and we shall 
call the vortex magnetic if cr^ > 0, electric for cr^ < 0, and chiral if o"^ = ^. 

d. It keeps its form under gauge transformations (12.121) generated by U=exp{ — |rr^} 
with constant F, whose effect is to rotate the field amplitudes, 

r r r r 

/i fi cos 2-/2 sm -, /2 ^ /2 COS 2 + /i 

ui ui COS r + M3 sin r, M3 — i> M3 cos V — ui sin F, 

171 ^ f 1 cos r + t>3 sinF, f3 ^ f3 cosF — t>i sinF, (4.8) 

whereas u,v rest invariant. These transformations can also be written in a compact form 
using the complex variables: 

(/i+«/2) e^^(/i+z/2), (Mi+zn3) ^ e-*^(ni+m3), (^^i+^^^s) ^ e-^^(wi+zt;3), (4.9) 

This symmetry can be fixed by requiring that 

fi(0) = 0. (4.10) 

e. The ansatz does not change if we multiply the amplitudes u,ui,U3 by a constant 
dividing at the same time cTq, by the same constant. To fix this symmetry, we impose the 
condition 

U3{0) = 1, (4.11) 

which is possible if ^3(0) 7^ 0. 

It should be said that the gauge (14. 4p is very convenient for calculations, since everything 
depends only on p, but it is not completely satisfactory, because, as we shall see below, the 
functions v, Vi, V3 do not vanish at p = and so the vector fields are not globally defined in 
this gauge. This problem can be cured by passing to another gauge (Eq. (l5.6p ). 
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A. Field equations 

With the parametrization fl4.4p the U(l) equations fl3.6al) . fl3.6cp reduce to 

^\^{u + us)f^ + 2uJJ, + {u-us)f^Y (4.12) 

p(^) = ^{(^ + ^^3)/i' + 2Vi/2 + (^-^3)/2'}, (4.13) 
where the prime denotes differentiation with respect to p. The Higgs equations fl3.6ep become 

kpnr = { Y [(- + «3)^ + + ^ + -3)^ + vi] + ^{fi + fi - 1)} A 



+ ( Y««i + ^^^^i)/2, (4.14) 



The Yang-Mills equations ( ]3.6bp .( l3.6dp reduce to 



-^{fM[)' = -1 {v,u, - v,u,) t;3 + ^ [«i(/i2 + fl) + 2^/1/2] , (4.16) 
11 

-(p«' )' = + (t;iU3 - t;3Ui) + ^ [(^3 + m)/i + (% " ^i)/|] , (4-17) 

p 2 

P (t") = ("^1^3 - ^3Wi) W3 + y [^i(/i + /I) + 21^/1/2] , (4.18) 

P {j^ = -Cy^ (^l«3 - V^U^) «1 + Y [(^3 + + (173 - V)f^ . (4.19) 

In addition, a careful inspection reveals that, although we have set to zero the radial and 

the second isotopic components of the Yang-Mills field, the Yang-Mills equation (13.6dp with 
a = 2 and k = p is not satisfied identically but gives the condition 

A = 0, (4.20) 

where 

A ^ a'iu.u', - u,u[) + 1 M - vsv[) - g\hr, - hf[). (4.21) 
Differentiating this quantity and using Eqs. fl4.12p - fl4.19p we discover the following relation 

A' + ^ = , (4.22) 
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such that 

A = ^ , (4.23) 

where C is an integration constant. This is a first integral for the above second order 
equations, so that replacing one of them by Eg. (14.231) would give a completely equivalent 
system. We thus see that Eq. (l4.20p plays a role of constraint restricting the value of C. 
The radial energy density for axially symmetric fields is pT^ = p{Si + S2), where 



= + 1« + n-) + -i-,(.,n3 - .3^1)^ 



+ I [{{u + u,)h + ni/2)' + ((n - U3)f2 + u^hf] I (4.24) 



and 

^2 = ^ [-^V^ + ^« + ^^3^) J + /f + /^^ (4.25) 
+ [((^ + ^3)/i + vj2f + {{v - vs)f2 + v^f^f] + ^ {fl + /| - 1)^ . 
Using the above equations it is not difficult to see that p£i is actually a total derivative, 

p£, = '"-^^i^j-,uu'+-^{u,u\ + u,u',)^ , (4.26) 
in agreement with Eq. (13.251) . Next, one has 



+ ( ^ - (^^ + ^^3)f - (^ - ^3)f - vjj2 - ^ P ) (4.27) 



with 



B2 = v[+g'phf2, 

B, = v', + g'P-Ul-f2l 

B4 = f[ + Ypl(^ + vs)h + vj2], 

B, = f!, + ^^[{v-vs)f2 + v,h]. (4.28) 

It follows that if P > 1 then pS2dp is bounded from below by the integral of the total 
derivative in the second line in (14.271) . The bound is achieved for (3 = 1 if Bi = . . . = 
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^5 = 0. Solutions of these Bogomol'nyi equations automatically fulfill the five second 
order equations fli?T3l) - flirT5l) . fH?T8l) . fH?T9l) with = 0. However, unless for g = 0, the 
Bogomol'nyi configurations are not very interesting, since, as was discussed above, they do 
not approach the vacuum at infinity and their energy is infinite. In what follows we shall 
rather focus on solutions of the second order equations fl4.12l) - fl4.20l) for generic values of (3 
and g. 



In order to construct global solutions of equations f l4.12l) - fl4.20l) in the interval p G [0, oo) 
we shall need their local solutions in the vicinity of the singular points, p = and p = oo. 
We shall be considering fields that are regular at p = and approach the vacuum for p ^ oo. 
Let us first consider local solutions at small p. 



Expressions fl4.24p . fl4.25p for the energy density can be used to derive the regularity 
conditions at p = 0. For regular fields the energy density at the axis must be bounded, and 
so the coefficients in front of the negative powers of p in fl4.24p . fl4.25l) should vanish at p = 0. 
This implies that at p = one should have 



V. BOUNDARY CONDITIONS 



A. Boundary conditions at the symmetry axis 



V = V 




(5.1) 



and also 



V1U3 



V3U1 = 0, 



{V + V3)fi 



+ 



Vlf2 = 



{V -V3)f2 + Vifi = 0. 



(5.2) 



Using Eq.f l4.10 l ) these conditions reduce to 



V3U1 = 0, 

iv + V3)fl=0, 
{V - Vz)f2 = 0. 



(5.3) 
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Let us now remember that the azimuthal components of the vector fields should vanish at 
the symmetry axis for the fields to be defined there. In the gauge (14.41) this would require 
that f (0) = ^3(0) = 0, but imposing such a condition would be much too restrictive. A 
more general possibility is to perform a gauge transformation that gives <y9-dependent phases 
to the scalar fields and shifts the azimuthal components of the vectors, and only after this 
to impose the regularity conditions for the vectors. Since the two Higgs field components 
should be single-valued in the new gauge, their phase factors should contain integers. We 
therefore apply to the ansatz (14. 4p the gauge transformation generated by 

U = exp|^(r7 + 7/;r3)|, (5.4) 

with 

rj = (2n — z/)(^ + cTqX", ip = ftp — GqX'^ , (5.5) 
where n, v are two integer winding numbers. This gives 

W = {u{p) + 1 + tIu^{p) + r^K(p) - 1]} (5.6) 

e*"'^fi(p) 

^^^n-v-v[p)-rlv,{p)^T^\v-V;[p)\]d^, $= ' , 

^i{n-i/)ip+iaax"' j!rp\ 

with = \Jt^\J~^ = r^cosip — r'^sinip. We can assume without loss of generality that 
n > 1 and we shall see below that I < ly < 2n. The dependence of the fields on will be 
convenient in what follows but not essential for the regularity issue. 

The azimuthal components of the vectors in (15. 6p will vanish at the axis provided that 

v{0) = 2n-u, 1)3(0) = u, (5.7) 

and also if vi{0) = 0, as required by Eq. fl4.1UI) . This gives a larger set of the allowed 
boundary values than in gauge (14. 4p . In addition, the regularity of the scalar fields requires 
that /i(0) = and, unless for z/ = n, that /2(0) = 0. The conditions (15. Sp now reduce to 

z/Mi(0)=0, n/i(0)=0, (n-z/)/2(0) = 0. (5.8) 

Summarizing, the boundary conditions at p = are given by 

Mi(0) = 0, %(0) = 1, t;i(0) = 0, vsi0) = i^, vi0)=2n-iy, /i(0) = 0, (5.9) 

while m(0) and /2(0) can be arbitrary if i/ = n, whereas for 7^ n one should have /2(0) = 0. 
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We can now work out the most general local series solutions of the equations for small p. 
The corresponding analysis is presented in the Appendix A, the result is 



u 


= ai + . . . , 


Ui 


= a2p'' + 


U3 


= 1 + ..., 


Vl 


= Oip^^'), 


V3 


= u + asp^ + 


V 


= 2n — u + p'^ + 


fl 


= 05 p" + . . . , 


f2 


= gp'"-"' + • • • , 



where 01,02,03,04,05 and q are six integration constants. The dots here stand for the 
subleading term. As explained in the Appendix, all such terms will be taken into account 
in our numerical scheme. 



B. Boundary conditions at infinity 

We want the fields for p ^ 00 to approach vacuum configurations with zero energy 
density. The energy density £^1 +£^2 is given by Eqs. (14. 241) . (14. 251) and can be decomposed into 
the kinetic energy part containing the derivatives of the field amplitudes and the potential 
energy part that contains no derivatives. The potential energy is a sum of perfect squares 
that will vanish if only each term in the sum vanishes. For example, the last term in (I4.25P 
will vanish if only + /| = 1 and so 

/i(oo) = cos |, /2(oo) = sin |. (5.11) 

The value of 7 is an essential parameter characterising the solutions, we shall call it vacuum 
angle. It is worth noting that the global symmetry (14. 8 p acts as 7 — 7 + F, but this cannot 
be used to change the value of 7, since this symmetry has already been fixed by the condition 
(l4.1Up at the origin. 

At the same time, just for discussing the local solutions at large p, it is convenient to 
temporarily set 7 = to simplify the analysis, since afterwords one can apply the transfor- 
mation (14. 8 p to restore the generic value of 7. We therefore assume for the time being that 
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at infinity one has 

fl = 1, f2 = 0. 

With this, the potential energy part of Eqs. (14.241) . (14.251) will vanish if only 

Ml = f 1 = 

and 

u = -U3, V = -V3. 

Under these conditions the field equations (I4.12p - (l4.19p reduce to 



ipu' 



0, 



^1 =0, 

P. 



(5.12) 



(5.13) 



(5.14) 



(5.15) 



from where 



u = ci + Qln{p), v = C2 + Ap'^, (5.16) 

where Ci, C2, Q, A are integration constants. One should set A = since otherwise the kinetic 
part in the energy density diverges at infinity. However, one should keep the logarithmic 
term in (I5.16p . since the energy density approaches zero for p — 00 even if Q 7^ 0. As a 
result, the conclusion is that for large p the fields approach the following exact solution of 
the field equations. 



;i-r=^)((ci + Qlnp)a„cix"-C2ci¥p), $ 



(5.17) 



Computing the A, Z fields defined by Eq. (12.231) gives 



Af,dx^ = ^ ((ci + Qlnp) aadx" - C2 dip) , = 0. 



(5.18) 



This is the electromagnetic Biot-Savart solution describing fields outside an uniformly 
charged electric wire. The charge and current are obtained with Eqs. (l3.13p . (l3.14l) . 



Q 



dT^ad'^x = - j> eksdkAadx' = -2n , 



99 



(5.19) 



assuming that the charge/current distributions are smooth inside the wire so that the Gauss 
theorem applies. In addition, there is also a magnetic flux along the wire. 



= i Akdx^ = ~—C2. 
J 99' 



(5.20) 
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In what follows we shall construct smooth vortex solutions in which the wire is represented 
by a regular distribution of massive non-linear fields, while in the far field zone the massive 
modes die out and everything reduces to the Biot-Savart configuration (15.171) . The energy 
density at large p will then be proportional to Q"^ / and the total energy per unit vortex 
length will be logarithmically divergent, just as for the ordinary infinitely long electric wire. 
The energy of finite vortex pieces, as for example vortex loops, will be finite. 

Let us now consider the deviations from the electromagnetic configuration (15.171) within 
the full system of equations (I4.12I) - (I4.19I) in order to determine how the solutions approach 
their asymptotic form. The corresponding analysis in the linear approximation is carried 
out in the Appendix B and the result is given by Eq.l lB.lSl) . We only need to apply to 
Eq.l lB.lSl) the phase rotation (14.81) to restore the generic value of 7 in Eq. (l5.11I) . This gives 
the large p behaviour of the solutions, 

M = ginp + 01 + ^6-™^"+... 
VP 

v = C2 + Cig'^^e-'^^P + ... 
Ml + ^3 = e-'^ [^e-^"^""" + 2 [-Qlnp - ci + ^ e"'"^^]) + . . . 

V, + iv^ = e-*^ [cs^pe-^^^'''' + z [-C2 + c.g'' ^p e-"^^^]] + ... 

+ = 1 1 + ^ e""^"" + i — e-^ '"-^A + • • • (5-21) 

[ Vp Vp J 

These local solutions contain 10 independent parameters Ci, . . . , Cg, 7, Q and they approach 
(modulo the 7-phases) the Abelian configuration (15.171) exponentially fast as p ^ 00, with 
the rates determined by the masses mz,mH and rria. Here 



rrio 



^ml + a\Q\np + ciy (5.22) 



reduces to the W boson mass = g/V^ if cr^ = 0. For cr^ 7^ this could be viewed as 
the W boson mass 'screened' (cr^ < 0) or 'dressed' (a^ > 0) by the interaction of charged W 
bosons with the long-ranged Biot-Savart field (I5.17p . We also note that the solutions (I5.2ip 
satisfy the field equations (I4.12p - (l4.19p but not the constraint (I4.20p . In fact, the latter has 
already been imposed on the local solutions (I5.10p at small p, and since it 'propagates', it 
will be automatically enforced by extending the local solution (I5.10p to large values of p. As 
a result, it would be redundant to impose it again. 
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We can now outline our strategy for solving the field equations. The local solutions for 
small and large p are given, respectively, by Eqs. fl5.10p and (15.211) . Within the numerical 



multiple shooting method [33j we extend these asymptotic solutions to the intermediate 
values of p and match them. We use the multi-zone version of the method {m], with many 
zones whose number and sizes are adjusted to make sure that for large p the asymptotic 
solution (15.211) is a good approximation. In order to illustrate the procedure, let us describe 
the simplest case where there are only two zones and one matching point. To match the 
solutions at this point we have to fulfill the 16 matching conditions for the 8 field amplitudes 
and for their first derivatives via adjusting the free parameters in Eqs. (l5.10i) . (l5.2ip . It is 
therefore essential to have enough of parameters, since otherwise the matching conditions 
will be incompatible. Counting the parameters we discover that there are 17 of them: 6 
integration constants ai, . . . , 05 and q in (15.101) . then 10 integration constants Ci, . . . , Cg, Q, 
7 in (15.211) . and finally o"^. This is enough to fulfill the 16 matching conditions and to have 
one parameter left free after the matching. This parameter will label the resulting global 
solutions in the interval p G [0, 00). It is convenient to choose it to be g - the coefficient in 
front of the amplitude /2 in (I5.10p . 

Before integrating the equations, some preliminary considerations are in order. 

VI. CONSERVED QUANTITIES 

Keeping only the leading terms, the boundary conditions for the field amplitudes for 

^ p ^ 00 are 

ai M — » ci + Q In p , 
2n — u *—v ^ 

^ ui ^ — (ci + Q In p) sin 7 , 

1 ^ -(ci + Qlnp) C0S7 , 
^ fi ^ — C2 sin 7 , 

z/ ^ — > — C2 cos 7 , 

7 

OsP" ^ /l ^ COS - , 

(6.1) 
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The knowledge of these boundary conditions is sufficient to calculate most of the global 
physical quantities associated to the solutions. 

A. Fluxes and current 

Since the fluxes and currents are gauge invariant, one can compute them in the gauge 
( 14. 4p . in which case the calculations are simpler. The non-zero components of the field 
tensors then read 

B pa CTq,?/ , B pip V , 

W^„ = aX, W^^ = -<, a = 1,3; 

W^v, = (^1^3 - U'iVi)aa- (6.2) 
Introducing the function f2(p) defined by 



o /l /I • o fa o\ 

cosn = .2 , f2 ' smfi = -2-— 72, (6.3) 



such that fi(cx3) = 7, one can define three orthogonal unit isovectors 

= (sin f2, 0, cos fi) , 
/c" = (cosfi, 0, - sinfi), 

r = (0,1,0). (6.4) 

Here corresponds to the definition fl2.27p . The Nambu electromagnetic and Z fields f l2.26p 
then read 

Fpa = { — u' — — {u[ sin + Mo cos n) ) CTq, , 

V^' 9 J 

g g' 

Fp^ = f ' H (f ^ sin J7 + v'^ cos fi) , 

9' 9 

Zpa = iu' + u\ sin + M3 cos VL)aa , 

Zp^p = —v' — v[ sin Q — v'^ cos Q , (6.5) 
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whereas the 't Hooft (12.281) fields have the total derivative structure, 

= ( ^ u - — (mi sin + M3 cos fi) ) a, 
\9 9 J 



a 1 

/ \ I 



^Z, = i - V + — (vi sinQ + V3 cos Q) 
\ 9 9 

^fa = (u + Ui sin + U3, cos fi)Va , 

= — (t; + f 1 sin + t>3 cos fi)' . (6.6) 

The electromagnetic current density is 

J" = d^F^"" = ^ {pFP"")' (6.7) 

and integrating over the p, if plane gives the same result for the both definitions of the fields, 

la = -27r / {pF^Jdp = a„ . (6.8) 

JO 99 

For most of the solutions considered below the vector 0"^ = (ctq, 0-3) = (ctq, is spacelike, so 
that there is the rest frame where cTq, = a 6^. The restframe value of the current is la = 
with X = —2TTQa/{gg'). The restframe components of the electromagnetic field strength 
and Z field strength are 

Bz ~ ^pf^ Ppzi Hz ~ ^pipi Zpz, (6.9) 

such that the magnetic fiuxes are 

/»00 /"OO 

^p = 2n pB^dp, ^^ = 27r/ pH^dp. (6.10) 
Jo Jo 

In the Nambu case these integrals have to be computed, while in the 't Hooft case they 

evaluate to 

271 

*f = -{c2 + {iy- 2n)g^ + vg'^ cosfi(O)), 

99 

*| = 27r(2n + z/(cos^](0) - 1)). (6.11) 

One can also define the W-condensate components by projecting the restframe values of 
W^j, onto the unit isovectors defined in (16.41) . 

wz = -^ k'^w;^, = -pw;„ = ~ r^%. (6.12) 
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B. Energy, momentum, angular momentum 

Using Eqs. fl4.25l) . fl4.26l) and the boundary conditions (16. II) . the total energy evaluates to 

0-2 , ^2 roc 

E = n ° Qu{oo) + 2n / p^a^p. (6.13) 
9^9^ Jo 

Here the first term on the right is the contribution of the total derivative in Eg. (14.261) . 
Comparing Eqs. (l3.2ip . (l3.25l) and (I4.26P one can see that the same total derivative determines 
the value of the momentum, 

P = 2n^Qu{oo). (6.14) 
9 9 

We finally notice that, since in the axially symmetric case there is an additional Killing 
vector -ft'(<^), associated to this there is the conserved angular momentum 

M = I TlK^^^d'x = I T^d'x. (6.15) 

This also has a total derivative structure [32], since one has 



= a° (^-^ u'v' + ^ + 2a°3fJ {^\u + UaT''){v + t;,r^)$) (6.16) 

with a = 1,3. Using the field equations fl412l) . fj416l) . fl4Tn) gives 

pT° = + ^ + ^3«^)) ' (6.17) 

so that 

M = 2.^^Qv{oo) = ''-^^. (6.18) 
9 9 9 9 

Summarizing, the fluxes, current, energy, momentum and angular momentum are given 

by the above expressions, and in most cases they are totally determined by the boundary 
conditions (16. ip . 

VII. KNOWN SOLUTIONS 

The only known solutions of Eqs. (14. 120 - 04.200 for generic values of q, g' are the embedded 
ANO vortices. These exist in two different versions, called Z strings Sj and W strings 
corresponding to two nonequivalent embeddings of U(l) to SU(2)xU(l). 
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A. Z strings 

These solutions are obtained by setting u^, = —u = 1 and f2 = Ui = vi = after which 
equations (l412l) - fH:20l) reduce to 

p[-p) =^-^{v + vs)fl (7.1b) 



-{Pfir =[-^2(^ + ^3)' + f (/i - 1) ) /i , (7.1a) 



With 



p(|) =^(-3 + ^^)/^. (7.1c) 



= ^ (fi'^^ - S-'^^a), ^^ANo = ^(^^ + f3), /ano = /i (7.2) 



these equations reduce to the ANO system 



kpfLo)' = + f (/a\o - 1)) /ano , (7.3a) 



?/ \ ' 1 

'^ANO 1 _ -"^ i-2 



P J 2 

and to 



O /ano "^ANO) (7.3b) 



Solutions of the ANO equations (17. 3p can be found numerically, they comprise a family 
labeled by n = 1,2,... with the boundary conditions 

/ano -^1, Wano (7.5) 

for ^ p ^ oo. The only bounded solution for A is the constant whose value is determined 
by the boundary conditions (16.11) at the origin, A = ng"^ — u/2. 
Summarizing, Z string profiles are given by 

u = —1, V = 2(7'^(t>ANo — n) + 2n — u, Ui = 0, u-^ = 1, 
vi = 0, V3 = 2g^{vANo -n) + jy, fi = /ano, /2 = 0, (7.6) 

from where one can read off the values of the asymptotic parameters (16. ip . 

ai = ci = —1, (5 = 0, C2 = 2ng^ — u, = f2o = 7 = 0. (7.7) 
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Calculating the fluxes gives the same result for the Nambu and 't Hooft definitions, 



(7.8) 



while the current, momentum and angular momentum vanish, so that the energy is finite. 
The dependence of the amplitudes in (17.61) on the second winding number v is pure gauge 
for these solutions, since in the gauge (15. 6p it disappears: 

/ 



e^"'^/AKo(p) 




(7.9) 



B. W strings 



These are obtained by setting in Eqs. (I4.12I) - (I4.20I) 



M = u = til = Ms = f3 = /2 = 0, vx = 1vp,^o{p), /i = /ano(p)- 



(7.10) 



Writing this solution in the gauge (14.41) and applying the gauge transformation generated 
by U= exp(m(y5r^) gives 



W = 2r\n - Mako(p)) d^, 0i + z02 = /ANo(p)e^"^ 



(7.11) 



which is the parametrization used in Ref. 



. However, throughout this paper we are using 



different gauge conventions, in particular we impose the condition (14.101) . We therefore 
return to (I7.10p and perform the global gauge rotation (14. 8 P with F = — 7r/2, which gives 

1 



u = v = Ux=U':^ = Vx = ^, W3 = 2t;ANo(p), /i = -/2 = /ano(p)- 
From here we can read off the values of the asymptotic parameters. 



(7.12) 



V = 2n, Q = Ci = C2 = 0, 'J = Q = Qq = 



IT 



(7.13) 



Calculating the fluxes (16.10p . (l6.1ip gives, in the both definitions, 

47m _ Attu g' 



(7.14) 



v^' V2 9 

These values are not the same as in the Z string case, in particular the electromagnetic flux 
is now different from zero. W strings are therefore physically different from Z strings. 
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VIII. SMALL CURRENT LIMIT BOUND STATES AROUND Z STRINGS 

Our goal is to construct solutions of equations fl4.12p - fl4.20p more general that the em- 
bedded ANO vortices. Our strategy was briefly summarized above: we numerically extend 
the asymptotic solutions (15.101) and (15.211) to the intermediate region and impose there the 
16 matching conditions, which can be fulfilled by adjusting the 17 free parameters in the 
local solutions. This leaves one extra parameter, g, which determines the value of the lower 
component of the Higgs field at the origin. 

The matching conditions are resolved iteratively, within the standard method described 



in 3^. A good choice of the initial configuration is then important, since otherwise the 
iterations will not converge. The idea is therefore to start the iterations in the vicinity of an 
already known solution, for which we can only choose the Z string (it is unclear whether W 
strings can also be used). Let us therefore choose the Z string solution (17. 6p as the starting 
point of our analysis. One has in this case /2(p) = so that q = 0. 

Suppose now that < g ^ 1. It is then natural to expect the solution to be a slightly 
deformed Z string. It can therefore be represented in the form 

U = —1 + Su, V = Vz + Sv, Ui = Sui, Us = 1 + SU3, 

vi = 6vi, V3 = vz3 + 5v3, f^ = f^ + 5f\, /2 = 5/2, (8.1) 

where = 2g''^{vj,^o - n) + 2n - u, v^s = 2g'^{v^^o - n) + u, = /ano and Su,...,5f2 
are small deformations. Inserting this to Eqs. fl4.12p - f]4.19p and linearising with respect to 
the deformations, the resulting equations split into three independent groups: two coupled 
equations for 6u, 6U3 plus three coupled equations for 6fi, 6v, 6V3 plus three coupled equations 
for 5/2, 6ui, 6vi. The first two groups do not admit interesting solutions, so that we can set 
6u = du^ = 6fi = 6v = 6V3 = 0. 

The remaining three equations describe a slightly deformed Z string: W + SW, $ + 5$, 
where W and $ are given by Eq. f]7.9p and 

r+ / \ 

SW = — e'^[6ui{p)(T^dx'' - 6vi{p)dyD] + h.c. , 5$ = , (8.2) 

2 \^ e*("^-^)5/2(p) J 

with ip = VLf — (TaX" and r+ = + ir"^. Here h.c. stands for Hermitean conjugation. The 
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in(l + p) 



Hi + p) 



FIG. 1: The bound state solutions of Eqs. (j8.3p with sin^ = 0.23, /3 = 2 for n = u = 1, 
0-2 = 0.5036 (left) and forn = 1, = 2, = -0.0876 (right). 



amplitudes 6/2, Svi, 6ui fulfil the equations 



7 ip^ny 



4p2 



+ ^ (/,^ - 1) + + if^M - ^/z5« 



1 

V,2, . 



cr 



2p2 



1 5 



P 



P 



2 /z + 0" ) - o- t;z3 + g Vzfz 5/2 



(8.3a) 
(8.3b) 
(8.3c) 



This can be viewed as a spectral problem with the eigenvalue cr^. In fact, by suitably 
redefining the variables one can rewrite the equations in the form 



•4) 



where \l/ is a 3 component vector and is a symmetric potential energy matrix depending 
on the parameters of the background Z string solution, f3, 6^, n, and also on u. Although the 
dependence on u for the Z string background is pure gauge, this is not so for the deformations 
of this background, since for different values of u in (18.31) one obtains different results. 

We are interested in localized, bound state solutions of Eqs. (l8.3p . These equations admit 
a global symmetry 



5/2 ^ 5f2 + - fz, 5ui 5ui + r, 5vi ^ 5vi + T 



(8.5) 
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0.5 1 1.5 2 0.5 1 1.5 2 



FIG. 2: The eigenvalue cr^ of the spectral problem (j8.3p against /? for n = u = 1 (left) and for 
n = 1, 1/ = 2 (right) for several values of g''^ = sin^ 9w 

which is merely the linearised version of (]4.8p . Using this symmetry one can impose the 
gauge condition 

5/2(00) =0, (8.6) 
and then the local behaviour of the solutions at large p can be read off from Eq.f l5.21|) . 

6u, = ^e-"''^p + ..., 5v, = Cs^e-"''^P + ..., 5/2 = 0,^6-"^^^ + .... (8.7) 
VP 

Here 

ml = ml + a\ (8.8) 

while the logarithm present in Eq. (15.221) does not appear in the linearized theory. The local 
behavior at small p can be obtained from Eq. fl5.10p . 

6ui = a^p" + r + . . . , 6vi = Tvs + ..., 6h = qp^''-'^ + ^ /z + . . . , (8.9) 

where we have included the free parameter T, since the symmetry (18. 5p is now fixed at 
large p and not at p = 0. In fact, the gauge condition (18.60 imposed at infinity generically 
implies that 5t'i(0) 7^ 0, which does not agree with the previously adopted in (I4.10p condition 
fi(0) = 0. However, the advantage of this gauge is that Z string deformations are clearly 
localized around the string core (see FigH]). One can always apply the symmetry (18.50 to 
return to the gauge where 5t'i(0) = 0, but then the amplitudes 5ui,5vi,5p will not vanish 
at infinity. 
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The local solutions fl8.7p . fl8.9p contain 6 free parameters: C5, Cg, Cg, 04, F, g. One of them 
can be absorbed by the overall normalization and so there remain only 5, but together with 
0"^ their number is again 6, which is just enough to fulfill the 6 matching conditions to 
construct global solutions of three second order equations (18.31) . The solutions obtained are 
characterized by the following properties. 
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FIG. 3: The eigenvalue of the spectral problem (j8.3p against sin^ for /? = 2. The values of n 
and V are shown as (n, v\ In the o"^ < region the curves terminate when the condition (jS.lOp is 
violated. The chiral solutions with o"^ = are possible only for special values of d^^. 



For given /3 and for small enough sin^^ one finds 2n different bound states labeled by 
z/ = 1, 2, . . . , 2r;, with the eigenvalue = <J^{P, 6*^, n, v) (see FigH]). These solutions can be 
interpreted as small deformations of Z strings by a current Jq, ~ cTq,, although the current 
itself appears only in the next order of perturbation theory, li (3 > 1 then n of these solutions 
always have o"^ > 0. For n — 1 solutions the eigenvalue cx^ changes sign when increases, 
and there is one solution with < for any 6^ > Q (see Figs l2][3|) . If is negative then it 
cannot be too large, since solutions are localized if only the mass (18. 8p is real, which requires 
that 

(8.10) 



2 ^ 2 
cr > -m^ . 



Every Z string admits therefore 'magnetic' (cr^ > 0) and, unless is close to 7r/2, 'electric' 
(o"^ < 0) linear deformations. The 'chiral' (cr^ = 0) deformations are not generic and possible 
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■ iH 

CO 




FIG. 4: The a^{P, 9^, n^v) — curves for several values of (n, v). Curves with v <- n are confined 
in the /3 < 1 region, while those with v > n extend to larger values of /?. 

only for special values of /3, such that 

a\(3,e^,n,u) = 0. (8.11) 

This condition determines a set of curves in the /?, 6'w plane (see FigJlj), let us call them 
chiral curves. They coincide with the curves delimiting the parameter regions of Z string 
stability [9]. This can be explained by the fact that a simple reinterpretation of the above 
considerations allows us to reproduce the results of the Z string stability analysis 9|. Let 

((5W, (5$) ~ exp{±z(aox° + asx^)} (8.12) 

be the Z string perturbation eigenvector (18. 2p for an eigenvalue o"^ = cr| — ctq. One has 

(8.13) 



(To 



from where it follows that if > then choosing a-^ < a gives exponentially growing in time 
solutions, that is unstable modes, since the frequency Uo is then imaginary, while choosing 
(73 > o" (as was done above) gives stationary Z string deformations, since (Tq is then real. If 
(T^ < then (To is always real. The sign of therefore determines if the negative modes are 
present or not. As a result, the (t^(/?, ^w, n, z/) = curves separate the parameter regions in 
the (3, 6'w plane in which the z/-th deformation of the n-th Z string is of the electric/magnetic 
type, and at the same time the regions where the z/-th Z string perturbation mode is of 
stable/ unstable type. 
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IX. GENERIC SUPERCONDUCTING VORTICES 

Having constructed the slightly deformed Z strings in the linear theory, we can promote 
them to solutions of the full system of equations f l4.12l) -( l4.19l) . first for small q. In the 
fully nonlinear theory all 8 field amplitudes deviate from their Z string values, and not only 
■^1; "^1) /2- However, if the deviations (5f i, 5/2 scale as g for g -C 1, those for the remaining 
5 amplitudes scale as q^. 

Starting from small values of g, we increase q iteratively, thereby obtaining fully non-linear 
superconducting vortices. Equations (14.12p -( l4.19p can be viewed in this case as a non-linear 
boundary value problem with the eigenvalue cx^. It turns out that only the magnetic (a^ > 0) 
and chiral (a^ = 0) solutions are compatible with the boundary conditions at infinity. The 
main difference with the linearized case is that, once all fields amplitudes are taken into 
account, the logarithmic term appears in the expression fl5.22p for the effective mass, 

ml = ml + a\Q\np + cif. (9.1) 

This implies that for < one has the 'tachyonic mass' < at large p, in which case 
the asymptotic solutions fl5.2ip oscillate and do not approach the vacuum at infinity (this 



was not realised in Ref. 23|]). The energy then diverges faster than In p and the interpretation 
of this is not so clear, since such solutions cannot be considered as field-theoretic models of 
electric wires or charge distributions. We shall therefore not consider the < solutions 
in what follows. 

For the magnetic solutions with > the problem does not arise, since the logarithmic 
term in the effective mass (19.11) only improves their localization. As a result, the magnetic 
solutions do generalize within the full, non-linear theory. The described above multiplet 
structure of the perturbative solutions persists also at the non-linear level, up to the fact 
that only solutions with cr^ > are now allowed. For n = 1 one finds one such solution, 
with z/ = 1, while for n = 2 there are already three: with u = 1,2 and, provided that 6^ is 
not too large, also with z/ = 3 (see Figj3]). The general rule seems to be such that for a given 
n there are 2n — 1 solutions, of which those with = 1, 2, . . . n exist for any 9^ while those 
with = n + 1, . . . , 2n — 1 exist if only is not too close to 7r/2, as shown in Figs J2]l3l 

The regions of their existence are delimited by the chiral curves in FigHl For each given 
curve one has o"^(/3, 6'w, n, z/) > in the region helow the curve. Let us call it allowed region, 
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it corresponds to magnetic solutions that generalize within the full non-linear theory. The 
region above the curve corresponds to the electric solutions with < 0, we therefore call 
it forbidden region. The chiral curves with v < n are all contained in the upper left corner 
of the diagram where < 1, while those for u > n extend to the region where /3 > 1. It 
follows that ii V <n then the allowed region corresponds to any 6^ ii (3 > 1, and to that 
is not too close to7r/2if/9<l. If u > n then the allowed region corresponds to the lower 
part of the diagram located below the {n, v) chiral curve. 

The chiral curves in Fig H] are obtained in the limit of vanishing current, while for finite 
currents they remain qualitatively similar but shift upward. For large currents they seem to 
approach the upper boundary of the diagram, such that the allowed regions increase. There- 
fore, without entering too much into details, one can simply say that the superconducting 
vortices exist for almost all values of (3, 9^. 
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FIG. 5: Profiles of the superconducting vortex solution for (3 = 2, g'^ = 0.23, n 



1 and 



a = 0.373. The 'electric' amplitudes (left) show at large p the logarithmic growth typical for the 
Biot-Savart field. The remaining amplitudes (right) are everywhere bounded. 



The typical solution is shown in FigJS] These solutions describe vortices carrying an 
electric current - smooth, globally regular field-theoretic analogues of electric wires. At 
small p they exhibit a completely regular core filled with massive fields creating the electric 
charge and current. At large p the massive fields decay and there remains only the Biot- 
Savart electromagnetic field. Integrating the charge and current densities over the vortex 
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cross section gives the vortex 'worldsheet current' ~ cTq whose components are the vortex 
electric charge (per unit length), Iq, and the total electric current through the vortex cross 
section, 1^. Since = o"| — > 0, the vector la is spacelike, which means that there is 
a CO moving reference frame where ao = /q = 0. However, the current I3 cannot be boosted 
away, so that this is an essential parameter, which suggests the term 'superconductivity'. 
In what follows we shall denote by X the value of J3 in the comoving frame. We shall also 
call a twist, since in the comoving frame one has = a, which determines the 2;-dependent 
relative phase of the two Higgs field components in the ansatz (15. 6p . 

The solutions exist for any (3 > and for belonging to the region below the (n, i/)-chiral 
curve in FigHl Solutions depend on j3, 6^ and also on g, n, v so that is in fact a function 
of five arguments, cr^ = (T^(/5, 6'w, n, z/, g). In addition, reconstructing the fields (15.61) . every 
solution of the differential equations determines actually a whole family of vortices with 
fixed cr^ = (T3 — (Tq but with different values of (ctojCTs) related to each other by Lorentz 
boosts. For given /?, 6'w the superconducting vortices therefore comprise a four parameter 
family that can be labeled by n, z/, ao and q. These parameters determine physical quantities 
associated to the vortex: the electromagnetic and Z fluxes p and \E'^, the charge Jq per 
unit vortex length and the current /a, as well as the momentum along the vortex, P, and 
the angular momentum M. 




FIG. 6: The restframe current I (left) and twist a (right) as functions of the parameter q for the 
solutions with /? = 2, sin^ = 0.23. 

The parameter q measures the deviation from the Z string limit. Increasing q starting 
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from zero increases X. It turns out, however, that the dependence 1{q) is not one-to-one, 
since for q not exceeding a certain maximal value q^{(3, 6^,71,1/) one finds two different 
solutions with different currents, as shown in FigJHl while for q > q^, there are no solutions at 
all. The solutions thus show two different branches that merge for g ^ but have different 
behaviour for g — > 0, when the lower branch reduces to the currentless Z string, while the 
current of the upper branch solutions grows seemingly without bounds. It is worth noting 
that this behaviour is drastically different from that found in the model of superconducting 
strings of Witten (see Fig 121] in the Appendix C), where the maximal value of q corresponds 
to the 'dressed' currentless string. 

A similar two-branch structure emerges also for other solution parameters traced against 
q, as for example for cr(g) (see FigJE]). This suggests that solutions should be labeled not by 
q but by a parameter that changes monotonously, as for example the current I or twist a 
(see FigsITKH]). 



b 




2 4 6 8 10 12 2 4 6 8 10 12 



I I 

FIG. 7: The twist a (left) and vacuum angle 7 (right) against the current I for the solutions 
with (3 = 2, sin^ = 0.23. One has a ^ and 7 ^ 7r/2 when X grows, while the product al^ 
approaches a non-zero value, so that a ^ X^^. 

The electromagnetic and Z fluxes calculated according to the Nambu definition fl6.10p 
are shown in Figl9]as functions of the current. We see that in the Z string limit I the 
Z flux reduces to Aim while the electromagnetic flux vanishes, as they should. As explained 
above, in the full non-Abelian theory there is no reason for the fluxes to be quantized, and 
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FIG. 8: The parameters ci (left) and 02/ {2n — v) (right) against I for the n = 1,2 solutions with 
/? = 2, sin^ 6'w = 0.23. One has C2 ^ {2v - n)g^ for I ^ 00. 

so they vary continuously with I. 
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FIG. 9: The Nambu electromagnetic flux ^f/t^ (left) and Z flux ^'^/(47rn) (right) against the 
current for the vortices with (3 = 2, sin^ = 0.23. 



It is interesting that the fluxes computed according to the 't Hooft definition (16. lip are 
completely different and do not reduce to the Z string values in the Z string limit g — » 0. 
This is the consequence of the fact that the function 

f! - /I 



cosfi 



f? + fl 



(9.2) 
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approaches its limiting value for g — ^ non-uniformly. Indeed, for g = one has /2 = and 
cosfi = 1 everywhere. If g ^ 1 then cosVt 1 almost everywhere, apart from the origin 
where one has /i ~ p" and /2 ~ gpl""'^' so that cosfi = — 1 at p = 0. Using Eqs. fl6.lll) then 
gives the Z flux 



H 



Ann, 



q = 0, 

q^O, u = l,...,2n-l. 



(9.3) 



The formula (16. lip for the electromagnetic flux contains the parameter C2, whose value is 
known only numerically, but for g it reduces to C2 = 2ng'^ — u, which gives 



=0, 

*f = 47rz/^, 
9 



g = 0, 

g — i> 0, u = 1, . . . ,2n — 1. 



(9.4) 



We see that both electromagnetic and Z fluxes computed according to the 't Hooft deflnition 
are discontinuous for g — ^ and do not reduce to the Z string values when the solutions 
approach Z strings in the zero current limit. No such problem arises if one uses the Nambu 
deflnition of the flelds. 
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FIG. 10: The 'Coulombian charge' parameter Q (left) and the magnetic energy (right) against I 
for the solutions with (3 = 2, Sin — 0.23. 



As the vortices support the long-ranged electromagnetic fleld, their energy per unit length 
diverges logarithmically. The total energy consists of the electric and magnetic contributions. 
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E = El + E2, and it is the electric term that diverges, 

2 2 

Ei = 7r^^^±^Qu{oo), (9.5) 

since u{p) ~ Qlnp for p — oo. The 'Coulombian charge' Q is shown against the current in 
FigHO] and against 6^ in FigJT6l The magnetic energy E2 is finite and reaches its minimal 
value for q = q^, as shown in FigJTOl This suggests that the value g = is in some sense 
'energetically preferable'. Since this value determines the 'bifurcation point' between the 
two solution branches, one can expect that the stability of the solutions may change at this 
point. 



X. LARGE CURRENT LIMIT 



The superconducting strings in the U(l)xU(l) model of Witten exhibit the current 

n n n 

quenching - there is an upper bound for their current pj, y, 15|. In the Weinberg-Salam 
theory we generically do not observe this phenomenon. Large I corresponds to small a, 
and we were able to descend down to a ~ 10~'^ (after this the numerics become somewhat 
involved) always finding larger and larger values of X, without any tendency for this increase 
to stop. 

It seems that the reason why the current can be unbounded in the Weinberg-Salam 
theory is related to the spin of the current carries. The existence of the critical current 
in superconductivity models usually follows from the fact that too large currents should 
produce strong magnetic fields which will quench the scalar condensate and destroy the 
superconductivity. This is what happens in Witten's model where the current is carried 
by scalars, and so quenching the condensate quenches the current. Specifically, the current 
defined by Eq. flC.9l) of the Appendix C vanishes for 02 — > 0. In the Weinberg-Salam theory, 
as we shall see below, the very strong magnetic field also quenches the Higgs field inside 
the vortex. However, this does not destroy superconductivity, since the current is carried 
by vectors and so the Higgs field is not the relevant order parameter. In fact, the current 
defined by Eqs.f l2.26i) .f l3.13l) does not vanish for $ — 0. For large currents the vector boson 
condensate can be described by the pure Yang-Mills field, whose scale invariance implies 
that the current can be as large as one wants. In the Witten model this does not work since 
the current carriers are scalars which break the scale invariance. 
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At the more technical level the current quenching can be related to the 'dressed' current- 
less solutions. If they exist, then changing the twist a makes the system interpolate between 
the 'bare' and 'dressed' strings, and since both of them are currentless, the current passes 
through a maximum in between. As we shall see in the next section, 'dressed' solutions 
exist in the electroweak theory only for special values of /5, 6'w, which means that current 
quenching is non-generic. For generic parameter values there are no 'dressed' solutions, and 
when one starts from the 'bare' Z string and decreases the twist, the current always grows. 

For large currents a simple approximation for the solutions can be found that agrees very 
well with the numerics and suggests that the current can be arbitrarily large. 

The large current limit is characterized by the presence of two different length scales. 
One of them does not depend on the current and is determined by the Z boson or Higgs 
boson masses, which are both of the same order of magnitude (for physical values of (3), 
i?z ~ ~ ~ 1. Another scale is related to the 'dressed' W boson mass fl5.22p . which 
increases with current and becomes large for large currents, mo- ~ crQ ~ X, so that the 
corresponding length is small ~ m^^ ~ X^^. One can therefore expect the solutions to 
show a small central region where p < Ra that accommodates the very heavy field modes 
of mass For R^ ^ p ^ Rz these 'supermassive' modes die out but the other massive 
modes remain, while the asymptotic region where R^ < p contains only the massless modes. 
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FIG. 11: Solutions profiles for /3 = 2, g'^ = 0.23, n = = 1 and a = 0.008. 



The numerical profiles of the large X solutions (see Fig JTTj) essentially confirm these 
expectations, but also reveal additional interesting features. For large currents the vacuum 



46 



angle 7 tends to 7r/2 (see Fig|7j) so that one has /i = /2 = 1/a/2 for p ^ 00. However, as is 
seen in FigJTTl one has /i ~ /2 not only for large p but everywhere. In particular, there is 
a region where /i ~ /2 ~ 0, which means that the gauge symmetry is completely restored 
to SU(2) xU(l). In fact, at the origin one still has /i = a^p and /2 = 5' (see Eq. fl5.10p ) but 
the parameters q, are extremely small. In addition, one has in this region vi ~ au^ ^ 0, 
V ^ 1, au ^ const., while ^3 and Ui change very quickly, and in particular ^3 vanishes very 
rapidly. The latter corresponds to freezing away of the very massive modes. 

One can wonder how the existence of the very massive modes can be compatible with the 
complete symmetry restoration. However, these modes arise not via the Higgs mechanism 
but rather due to the screening by the current, which is somewhat similar to the colour 
screening in pure Yang-Mills systems. When V3 vanishes, the Higgs field starts to change 
and approaches its vacuum value, after which all massive degrees of freedom die out and the 
solutions reduce to the U(l) electric wire fl5.17p . 

It turns out to be possible to find a relatively simple approximation for the large current 
solutions that explains all their empirically observed features. Specifically, since the vacuum 
angle 7 approaches tt/2 for large currents, let us set 7 = 7r/2. Let us also restrict ourselves 
to the simplest case where n = u = 1. The boundary conditions (16. ip then become 

ai 4— -u — > ci + Q In p , ^ mi ^ — (ci + Q In p) , 1 ^ ^3 — > , 

l^f-^C2, ^ Vi ^ -C2 , 1 ^ f 3 -> , 

a,p^fi^l/V2, g^/2^i/v/2. (10.6) 



A. W-condensate region — SU(2) Yang-Mills string 

Let us approximate the solution at small p by 

fi = f2 = cru3 = vi = 0, au = aai, v = 1, aui{p) = \Ui{Xp), v3 = Vs{\p), (10.7) 
where A is a scale parameter. Inserting this to (I4.12p - (l4.20p the equations reduce to 



1 V'^ 



-{xU[y = ^Ui, (10.8a) 



X 

xC^^ =U^Vs, (10.8b) 
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with X = Xp. These equations admit a solution with the following boundary conditions for 

^ a; ^ oo, 

x + ... ^Ui{x) ^0.85 + 0.91 ln(x) + ... (10.9a) 
1 - 0.45x' + . . . ^^3(2;) ^ 0.32v^e°-°^"x-°-^^" + . . . (10.9b) 

for which V3 rapidly approaches zero (see Fig JT^ . Inserting Ui{x),V3{x) to fll0.7p gives a 
family of solutions distinguished from each other by values of the scale parameter A. This 
scaling symmetry arises due to the fact that, since the Higgs field is zero, the system is pure 
Yang-Mills. In fact, Eqs. (110.81) coincide with the Yang-Mills equations for the pure SU(2) 
system, 

Csui2) = -\'^%W'^'"' (10.10) 
provided that the Yang-Mills field is cylindrically-symmetric 

T^W^dx" = T^XUi{Xp)dz + T^Vs{Xp)dip. (10.11) 

We shall therefore call the solution shown in FigJT2] Yang-Mills string. 

The Yang-Mills strings are made of the pure gauge field and describe the charged W- 
condensate trapped in the region where the amplitude V3 is different from zero. Outside 
this region there remains the long-range Biot-Savart field generated by the condensate and 
represented by Ui. The right-hand sides of Eqs. (110.81) can be viewed as conserved (although 
not gauge invariant) current densities. Integrating over the p, (f plane gives the current along 
the z-direction, X ~ A, localized in the region of the size ~ 1/A where V3 is different from 
zero. 

Yang-Mills strings approximate the central current-carrying core of the electroweak vor- 
tex. Their scale invariance implies that there is no upper bound for the current X ~ A, since 
the scale parameter A can assume any value. 

B. External U(l)xU(l) region 

Let us now assume that the Yang-Mills string approximation is valid only in the central 
vortex core, for x < xq, with xq chosen such that Vs^xq) is sufficiently close to zero, for 
example Xq = 4 (see Fig JT2|) . This determines the upper value po = Xq/X. As we shall see, 
for large A the precise knowledge of Xq is unimportant. 
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FIG. 12: Left: solution of Eqs. (|10.8p with the boundary condition (jl0.9p and the function sinfi 
defined by Eqs. (f6:3l) . (jl0.26p . Right: the typical solution of Eqs. (jl0.17!l . 

For p > po we set 

vs = aus = 0, fi = f2 = f/V2, (10.12) 

while the remaining field amplitudes can be combined to define the electromagnetic and Z 
field variables, 



v = -{v+v,: 



Ua = cr{g\ - g'^m), Va = - g%,. (10.13) 

More precisely, using formulas of Sec. flVI Ap and Eq. fll0.12p . the electromagnetic and Z fields 
potentials are 

A, = j^B,-^Wl, Z, = B, + Wl (10.14) 

the Higgs field components are 0i = 02 = while the W-condensate components vanish. 
The original SU(2)xU(l) theory iKW reduces then to the U(l)xU(l) model 

CuiD.uw = --^F.uF^'' - \z,.Z^' + (2^,0)*P> - ^(101' - 1)' (10.15) 

with P^0 = [d^ — |Z^)0. This contains the free Maxwell theory and the Abelian Higgs 
model. The fields are gg' A^dx^ = UAdz + VAdip and 2Z^dx^ = Udz + Vd(f while = /. 
Eqs.f l4.12l) - fl4.20p then reduce to the Maxwell equations 

ipU'Ay = 0, =0, (10.16) 
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and to the equations of the Abehan Higgs model 



^(p/')'=(f/' + ^ + f(f -1))/, (10.17a) 

P{j)-lfy^ (10.17b) 

kpU'y = lfU. (10.17c) 
P 2 



The solution of the Maxwell equations (110. 16p is 



Ua = A + B\np, Va = C, (10.18) 

with constant A,B,C. Equations fll0.17p admit a one-parameter family of solutions with 
the boundary conditions 

u{po)^u^o, l^y^o, o^/^i (10.19) 

for pq <— p ^ oo (see Fig JT2l) . The numerics show that if U'{po) is large and positive then 
U{po) is large and negative and the derivatives /'(po) and V^'(po) tend to zero such that 
there is a neighbourhood of po where 

V^l/2, f^O, U^a + b\np (10.20) 

with the parameters a, b related to each other, a = a{b). As a result, 

f/(po) = a(6) + Mnpo, t/'(po) = b/po. (10.21) 



C. Matching the solutions 

The solutions obtained for p < po and for p > po should agree at p = po = Xq/A where 
the functions and their first derivatives should match. Matching u, ui and their derivatives 
gives four conditions that determine the free parameters ai, b, A, B in the above solutions, 

B = -0.9V^A, 
b = 0.45A, 
aai = 2a{b) - 0.85A - 26 In A, 
A = 51nA + /ai -O.SS^'^A. (10.22) 
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It is worth noting that Xq has dropped from these relations. Matching similarly v,Vi, f gives 
only one condition, C = g^, but the derivatives v',v[,f' do not match precisely, since they 
vanish for p < po but not for p > po- However, this discrepancy tends to zero when A grows, 
since U'{po) ~ then increases thus rendering the approximation (110.201) better and better. 
Similarly, although there is a discrepancy in values of V3, since it vanishes for p > po but 
not for p < Po, increasing A reduces this discrepancy as well. 

As a result, by matching the two solutions we can approximate the global solution in 
the interval < p < 00. The approximate solution depends on A and approaches the 
true solution when A increases. For A ^ 00 the discrepancy of values of V3,v',v[,f' at 
the matching point vanishes and so the approximate solution becomes exact. There only 
remains to relate A to the current. Comparing with (110. 6p we see that aQ = B (also aci = A 
and, as noticed in FiglHl C2 = g"^) and since the current X = —2'KaQ/{gg') it follows that 



A 



0.17-1. 
9' 



;i0.23) 



The approximate solution shown in FigHS] clearly exhibit the same structure as the true one 
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FIG. 13: The approximate solution for /3 = 2, g'"^ = 0.23, n = 1/ = 1 and A = 7.6, I = 24.43. 



in Figim Although it corresponds to not a very large value of the scale parameter, A = 7.6, 
its current is twice as large as compared to that in FigJTTl The true solution in FigJTTl 
shows a slowly growing amplitude U3 whereas in FigJlS] one has au^ = 0. This difference 
is a subleading effect, since in Fig{TT]the vacuum angle is not exactly 7r/2 and so there is 
a growing mode u-^ ~ cos 7 In p. For I ^ 00 one has 7 ^ 7r/2 and this mode disappears. 
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so that U3 becomes everywhere bounded, and multiplying it by a ^ will give au^ = 0, as 
shown in FigJT3l With some more analysis such subleading effects can be taken into account 
to determine, for example, 7(X), cr(X), g(X). 



D. Symmetry restoration inside the vortex 

Comparing FigsHT] and [13] reveals that the zero Higgs field region expands when the 
current increases. The inner structure of the large current vortex can be schematically 
represented by FigHH The vortex shows a large symmetric phase region of size ~ X where 
the Higgs field is very close to zero. In the centre of this region there is a compact core 
of size ~ X^^ containing the charged W-condensate and approximated by the Yang-Mills 
string. Outside the core there live the electromagnetic and Z fields. The symmetric phase 
is surrounded by a 'crust' of thickness ~ Rz where the Higgs field approaches the constant 
vacuum value while the Z field becomes massive and dies away. The region outside the 
'crust' is described by the Maxwell electrodynamics and is dominated by the Biot-Savart 
magnetic field generated by the current in the core. 




FIG. 14: The schematic structure of the iarge current vortex. 



This picture reminds somewhat the 'W-dressed superconducting string' scenario of Am- 
bjorn and Olesen [l^, briefly described in Sec lIIl Dl above. In this scenario Witten's su- 
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perconducting string produced at the GUT energies acts as the external source of magnetic 
field that changes the structure of the electroweak vacuum. There is a cylindrical layer 
around the string where the field falls within the limits Eq. (13.331) and where the Higgs field 
is non-trivial. Inside this layer the field is supercritical and the Higgs field is zero, while 
outside the field is undercritical and the Higgs field is in vacuum. All this is very similar to 
what is shown in Fig{TH This suggests that we have a self-consistent electroweak realization 
of the Ambjorn-Olesen scenario in which the Witten string is replaced by the Yang-Mills 
string of completely electroweak origin. The following estimates confirm this. 

First of all, the radius of the layer where the Higgs field is non-trivial should scale as 
X, and we shall now see that the crust radius scales precisely in this way. Indeed, inside 
the crust / is very small and one can neglect the right hand sides in f ll0.17cl) to obtain 
U = b\ia{p/ p^), where is an integration constant and b = 0.45A. As long as p ^ p^, the 
term in fll0.17ap is large and makes / stay very close to zero, since / can increase and pass 
through the inflection point if only f/^ is small. The latter follows from the fact that the 
second derivative of / with respect to In p vanishes if only the right-hand side of fll0.17ap 
vanishes, which is only possible if f/^ is small enough to be canceled by the negative term 
|(/^ — 1). Therefore, / can vary if only U ~ but as soon as / deviates from zero, it 
acts as the mass term for U to drive it to zero. As a result, / interpolates between zero 
and one close to the place where U vanishes, for p ~ p*, in an interval of size 6p ~ R^- 
Multiplying ( I10.17cp by p and integrating from po to inflnity gives 2b = — ppUdp and 
since the integrand is almost everywhere zero apart from a small vicinity of p^, it follows 
that b ~ v^P*- This gives the size of the symmetric phase region (the crust radius) 

p, = 0.28 X ^ — , (10.24) 
9' 

where the coefficient is found numerically. As a result, one can approximate the solution of 
Eqs.f fTTUTl) by 

/ ^ e(p - p.), f/ ^ 6 (1 - /) ln(p/pj, ^ 1(1 - /)(1 - (p/p,)2) (10.25) 

where 6(p — p*) is the step-function smoothed over a finite interval 5p ~ i?z. 

Let us now consider the distribution of the fields inside the vortex. The fields are 
defined by Eqs. (l6.5p . (l6.12p where there is a hitherto undetermined function Vt. Using 
Eqs. (l6.3p . (ll0.12p one finds that VL = tt/2 in the external region, while in the core fl is 
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not yet defined. Let us calculate the first order correction to the core configuration fll0.7p 
assuming that /i, /2 are small functions of x = Ap and that /i ~ a; and /2 = 0(1) for a; — * 0. 
Linearising Eqs. fl4.14p . fl4.15l) gives for large A, with /± = /i ± /2, 

iwiy.(4„)^,i±S);,,^;,. ,10.26, 

where p = —aai/{2X) ~ In A and f/i, V3 are given by Eqs. fllO.Sp . Integrating these equations 
and using Eq. fl6.3p determines VL inside the core and shows that it approaches 7r/2 very 
rapidly (see FigJT2l). 



140 I i i 1 1 1 1 i 1 2.5 




ln(l + p) ln(l + p) 

FIG. 15: The field distributions inside (left) and outside (right) the current-carrying condensate 
core of the vortex, for the same solution as in Figll31 



We have now all the ingredients to compute the fields fl6.5p .f !6rT2l) and the current density 
(16. 7p . The result presented in FigJlS] shows that the current, W-condensate, and the z- 
component of the magnetic field are trapped in the core, where the energy density is maximal. 
Outside the core there extend only the Biot-Savart magnetic field and the Z field, which are 
both massless as long as the Higgs field stays close to zero. At a distance ~ X away from 
the core the system enters the 'crust' region where the Higgs field approaches the vacuum 
value. This gives rise to a local maximum in the energy density and drives the Z field very 
rapidly to zero. Outside the 'crust' there remains only the Biot-Savart magnetic field. 

The most interesting feature is that the 'crust' region where the Higgs field is non-trivial 
exists only under certain values of the magnetic field. Outside this region, where the field is 
either too weak or too strong, the Higgs field is either in vacuum or vanishes, in agreement 
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with the predictions of [15]. At the same time, the magnetic field inside the crust does not 



quite fall within the limits (m^/e, m|/e) predicted in 15| for the homogeneous magnetic 
field with (3 = 1. The electroweak condensate in the latter case contains the W, Z and 
Higgs fields [15] • In our case the crust region does not contain the W field, since it develops 
the very large effective mass rrio- due to the interaction with the inhomogeneous Biot-Savart 
field, so that it can be different form zero only in the central core region. 

We also notice in Fig{T5]that, since Wp, do not vanish at p = 0, the W-condensate field 
strength is actually singular at the axis. However, this singularity is mild, since the energy 
density and the field potentials are globally regular. This singularity is absent for solutions 
with V > 1. 



XI. SPECIAL PARAMETER LIMITS 

The described above properties of the solutions remain qualitatively the same for the 
generic parameter values. However, for special values of /3, 9^ solutions can exhibit new 
features. Even though such values do not always belong to the physical region, where 
1.5 < /3 < 3.5 and sin^6'w = 0.23, it is instructive to consider them, since this helps to 
understand the structure of the solution space. 

For example, the limits 6'w — > or 6'w — ^ vr/2 are special, because either the hypercharge 
field or isospin field W"j then become massless and decouple. The long-range mode 
Q In p then disappears since Q vanishes (see FigJT6ll and the total energy of the solutions 



becomes finite. Solutions for 6^ = tt/2 have been studied in the literature [2l|, but it is 
worth reviewing them in the context of our general discussion. 

Another interesting case corresponds to the chiral solutions, which exist only for special 
values of 6'w given by a certain function 6'w(/3, n, u, q). This function can be determined by 
requiring the solutions to respect the o"^ = condition. Finally, we shall consider infinite 
Higgs mass limit, /3 — > cxd, in which case the theory reduces to the gauged CP^ sigma-model. 
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FIG. 16: Left: the 'Coulombian charge' parameter Q against 6^ for the solutions with (5 = 2, 
n = V = 1 and q = 0.2. Right: the twist a against q for solutions with n = u = \, j3 = 2. For 
g ^ the curves a{q) start at g = at their maximum, then go down towards larger values of q, 
but then always bend back to g = 0. If = then a decreases monotonously to zero at g = 1. 

A. Semilocal limit, = 7r/2 

Since the right hand side of the Yang-Mills equations (12.161) is proportional to g"^, one 
has ~ g"^ for small g so that the field and its kinetic term in the Lagrangian vanish 
for (yf — i> and the theory (12.101) reduces to 

C = -] B^^B^^- + (D^<l>)tD^$ - f (<|.t$ - 1)' , (11.27) 

with -D^$ = {d^ — I -B/i)$. This reduced theory is sometimes called semilocal, since it has 
the local U(l) invariance, while the SU(2) symmetry of the original theory now becomes 
global [2^. The perturbative mass spectrum of the model (13.41) can be obtained by setting 
(7 = in Eq. (l2.24i ): it contains a massive vector boson and a Higgs boson, as well as a pair 
of massless Goldstone scalars. However, the latter acquire a finite mass via the screening 
mechanism similar to that described above for the W bosons, so that all fields are massive. 

The field equations can be obtained by setting in Eqs. (l4.12l) - (I4.20I) the SU(2) field am- 
plitudes to be constant and equal to their values at the origin, ui = t>i = 0, U3 = 1, ^3 = z/. 
This solves the Yang-Mills equations (I4.16p -( l4.20p and implies the following values of the 
asymptotic parameters in (16. ip : Q = 7 = 0, ci = —1, C2 = —v. In the gauge (15. 6p the 
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Yang-Mills field then vanishes and the fields read 

e'"'^ fi (p) 

W = {u{p) + 1) ao.dx'' + {2n - ly - v{p))d^, $=| . |. (11.28) 



/2(P) 



Eqs.( H:T6|) - fH:20|) reduce to 



-^{pu'y =l(^{u + l)f! + {u- l)f^^ , (11.29) 

= K ^ " ''^^') ' ^^^'^^^ 

lipfiY = (y(« + 1)' + ^,iv + + ^(A^ + fl - 1)) (11.31) 



^(p/2)' = ^y(« - 1)' + ^(^ - + ^(/i + /2 - 1) ) /2- (11.32) 
The boundary conditions are obtained from Eqs. fl5.10p . (l5.2ip . 

ai + ... ^M-^ -l + ^e-'"2^ + ... . 



2n - z/ + a4P^ + . . . ^t; ^ -z/ + 045-'^ v^e"™^^ + . . . 

Vp 



asP" + . . . ^ /i ^ 1 + ^ e-"'^P 



gpl"-'^l + ...^/2^^e-"^-'' + ... , (11.33) 

where mz, mn are the same as in Eq. fl2.24p . while Eq. fl5.22p now gives 111^ = cr. Since all the 
fields are massive and approach zero exponentially fast, the vortex energy is finite. 

The simplest solution of Eqs. (11 1.29l) -f ll 1.321) is obtained for g = 0, this is the Z string 



(17. 6p restricted to 6'w = vr/2, in which case it is called semilocal string 201], 35|. Choosing 
g 7^ gives current-carrying vortices called twisted semilocal strings [21!], a typical solution 
being shown in FigJTTl These solutions exist only for f3 > 1 and for u = 1, . . .n. This is due 
to the fact that for 6^ = 71/ 2 the allowed region in Fig H] is (3 > 1 if u < n while for v > n 
the allowed region is empty. For given /3, n, v the solutions comprise a family that can be 
labeled by g G [0, 1). 

The current in this case is not local but global, related to the SU(2) invariance of the 
model (13. 4p . It is determined by in Eq. (l2.16p . 

//»CXD 
Jld'x = na^j^ {(u+l)f^ + il~u)f^}pdp, (11.34) 
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FIG. 17: Features of the n = v = l^(3 = 2 solutions in the semilocal hmit, g = {). The solution 
profiles for a = 0.1 (left) and the total energy E/{2TTn) against the restframe current (right). 

while the similar integrals of and vanish. The energy E = 2it /q°°(^^i + S2)pdp is finite 
and can be obtained form Eqs. fl4.24l) . fl4.25p by setting the SU(2) amplitudes to the specified 
above values and omitting terms with their derivatives. For (3 > 1 and g = there is a 
non-trivial lower bound for the energy determined by the total derivative in the second line 
of Eq.g27D, 

E > E^^ = 2n + ^)fi - - ^)/2}'^^P = 2vrn. (11.35) 

Jo ^ 

It turns out that the energy of the 'twisted strings' monotonously decreases with current (see 
Fig lTTI) and approaches the lower bound 27m when the current tends to infinity 2l|]. The 
twist (J also decreases with current, as for solutions with 7^ (see FigJT]). The parameter 
q increases monotonously and approaches unit value when the current tends to infinity 2l| , 
while for solutions with (7 7^ it always passes through a maximum and then tends to zero 
for large currents, as shown in FiglHl If one plots q against a (or against the current) for 
different then it turns out that the curves approach the limiting curve for (7 = pointwise 
but not uniformly, since the boundary condition at a = are different, depending on whether 
(7 = or (7 7^ (see FigIT6l). 

The lower energy bound -Emin in Eq. flll.35p is attained for /3 = 1 if the first order Bogo- 
mol'nyi equations Es = Q with Bs defined by Eq. (l4.28p are fulfilled. For g = Q there are only 
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three such equations, 

v'+^^M + f^)-l] = 0, 

+ ^ + i^)fi = 0, 

f2 + ^iv- u)h = 0. (11.36) 

Their solutions comprise a family labeled by g G [0, 1) 36j, they also fulfill the three second 
order equations flll.30p - flll.32l) with cr^ = 0, while the amplitude u can be obtained from 
Eq. flll.29p . Reconstructing the fields (111.280 . solutions with (Xa = are sometimes called 
'skyrmions', while those with ctq = =to"3 7^ can be called 'spinning skyrmions' [37!]. Their 



profiles are qualitatively similar to those for the twisted strings, but since rricr = cr = in this 
case, they show polynomial tails at large p. It turns out that the twisted strings with (3 > 1 
approach the 'spinning skyrmion' configurations for X ^ oo (modulo an infinite rescaling) 



B. Isospint limit, 0w = 

If g' = sin^w is small then the hypercharge field i?^ scales as g'^ so that for g' ^ the 
theory fl2.10p reduces to 

£ = -Iw^^^W""" + (D^<I>)^D'^$ - ^ ($^<l> - 1)^ , (11.37) 

with Df,<^ = (5^ - 1 r"W;j) $ . The SU(2) is now local while the U(l) is global. The pertur- 
bative mass spectrum contains three massive vector bosons with the mass rriz and a Higgs 
boson with the mass mn. The field equations in the axially symmetric case are obtained by 
setting in fl4.14l) - fl4.20l) u{p) = const., v{p) = 2n — u, which implies that the U(l) part of 
the field (15. 6p vanishes and that one has in Eqs. fl5.10p . fl5.2ip . fl6.ip ai = Ci = u, C2 = 2n — u 
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and Q = C3 = C4 = 0. The boundary conditions for the remaining field amphtudes read 

<— Ml — > —u sin 7 , 

1 ^ ^3 ^ —u cos 7 , 

<— f 1 — > (z/ — 2n) sin 7 , 

u ^ {u — 2n) cos 7 , 

7 

flsp" ^ /i ^ cos - , 
gpl"-^l-/2^sin|, (11.38) 

The simplest solution is obtained for q = 0, this is the Z string (17.61) restricted to 6^^ = 0. 
For q ^ one obtains current-carrying solutions. The typical solution is shown in FigJTSl 
These solutions relate to the lower boundary of the chiral diagram in FigJH which implies 
that they exist for any (3 > and for u = 1, . . . 2ri — 1. Their current is global, related to the 
global U(l) invariance of the theory (111.371) . it can be expressed in terms of J° in Eg. (12. 151) . 

//•oo 
J° d'x = naa {{us + u)fl + 2^1/1/2 - (^3 - n)/|} pdp. (11.39) 

As for all solutions described above, there are no apparent restrictions on its values. The 
energy density is defined by Eqs. (l4.24l) . (l4.25p with u,v restricted to the values specified 
above. The total energy is finite, it always increases with current (see Fig JTSj) . 

These = solutions can also be obtained from the generic electroweak vortices in the 
limit 6'w 0. Similarly the = 7r/2 solutions can be obtained in the limit Tijl. 
The generic field configurations then approach the limiting ones pointwise, although non- 
uniformly at large p, since they support a long-range field ~ Q In p, while in the limit one 
has (5 = 0. Since the electromagnetic current vanishes in the limit, one can wonder how the 
global currents (111.341) and (111.391) can be obtained via the limiting procedure. The answer 
is that they can be obtained as limits of the rescaled electromagnetic current Ia/{gg')- The 
following relations take place, 

g 

= ^ d^B,^ - - \ {D^n^)W^^ . (11.40) 

g g 
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FIG. 18: Features of the n = i/ = l,/? = 2 solutions in the isospin hmit, g = 1. The profiles for 
a = 0.5 (left) and the energy against the restframe current (right). The insertion shows that E/n 
increases with X faster for n = 1 than for n = 2, so that the vortices repel each other for small 
currents and attract otherwise. 

Here we have used Eqs. fl2J[5|) . fl2:T6D and the identity d^in^Wf^J = (D^n")iy;„ + n'^D^Wl^ 
where Do-n" = d^n"' + eafecW^n"^. Integrating the first line of these relations over the vortex 
cross section gives la- Integral of the second line reduces in the limit g' —>■ to the global 
current ( 111.39^ . since the total derivative term in this line gives no contribution, because 
the field W^^ becomes short-ranged in this limit. 

Let us now notice that, since the whole expression in (lll.40p has the total derivative 
structure, its integral will not change upon replacing the vector = ($V"$)/ ($^$) by any 
other unit vector with the same value at infinity. Let us therefore replace n°' by n°' = 6^ in 
the third line in (111.401) . One will have D^n"- = eabsW^^, which scales as for small g, so 
that the last term in the third line vanishes in the limit g ^ 0. The second term gives upon 
integration the global current (lll.34p (up to the sign), while the first term gives zero, since 
the field becomes massive in the limit. We therefore conclude that 

I^^'^ ^^rJ^^ -^f (11-41) 

for g' ^ and g ^ 0, respectively, so that the global currents can indeed be obtained from 
the local electromagnetic current via the limiting procedure. 
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C. Special chiral solutions 

Superconducting strings in the Witten model are usually constructed starting from the 
'dressed' cr^ = chiral/ current less string, viewed as the 'bare' ANO vortex stabilized by the 
condensate in its core. Increasing o"^ produces the current, while the parameter q decreases 
and for g ^ the solutions reduce to the 'bare' ANO vortex (see Appendix C). 

Superconducting vortices in the Weinberg-Salam theory are obtained in the opposite way, 
starting from the 'bare' ANO vortex (Z string) and then decreasing the twist a. One does 
not generically find chiral solutions in this case, since the limit cr^ ^ corresponds to infinite 
and not zero restframe current. 

However, chiral solutions can be obtained for special values of the parameters. As was 
discussed in Sec JVIIIt for g <^ 1 they exist if /5, 6^, n, v belong to the chiral curves shown in 
FigHl These curves are obtained by solving the linear eigenvalue problem (18. 3p . If q is not 
small, then one should solve the full system of equations (I4.12p - (l4.20p under the condition 

cr2(A^w,n,z/,g) =0. (11.42) 

This condition reduces by one the number of free parameters in the local solutions (IS.lOp . 
(I5.2ip . so that in order to perform the matching one has to consider as one of the shooting 
parameters. For given n, q this determines curves in the (/5, 6'w) plane. For g ^ these 
curves reduce to those shown in FigHJ while for g 7^ they look qualitatively similar but 
shift upwards. As a result, given a point j3, 9^ in an upper vicinity of a curve in Fig H] one can 
adjust g such that the shifted curve will pass through this point. This fine tuning determines 
the values of /3, 6'w, g for which there is a solution of Eqs. (l412il -( !419ll with = 0. The 
chiral solutions are therefore not generic but exist only for values of /3, 6'w which are close to 
the chiral curves in FigJH 

Reconstructing the fields (15.60 . a given solution with cr^ = cr| — ctq = determines a 
family of chiral vortices labeled by (Tq = ±o"3, different members of this family being related 
by Lorentz boosts. These vortices carry a chiral current and support a long-range field so 
that their energy is infinite. However, there is a distinguished solution with o"o = as = 
for which the current is zero and the energy is finite. This special chiral solution is not a Z 
string or its gauge copy, since it has ^ a 7^ (see FigJTOl). 



Such solutions remind of the 'W-dressed Z strings' discussed some time ago [lO|, lll |. 



By analogy with the 'dressed' vortices in the Witten model, these were supposed to be Z 
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FIG. 19: Profiles of the special chiral solution (left) and Z string (right) for n = 4, i/ = 7, /3 = 2, 
= 0.23. 

strings stabilized by the W-condensate in the core. However, it seems that this stabilization 
mechanism does not work in the electroweak case because Z strings are non-topological and 
can unwind into vacuum [3]. The numerical search for the 'W-dressed Z strings' gives a 



negative result [ij] and so it seems they do not exist. Fortunately, as we have seen, the 
'W-dressed Z strings' are not necessary for constructing the superconducting vortices. 

The special chiral solutions cannot be considered as the 'dressed' Z strings, first because 
they do not exist for generic values of the parameters, and secondly because they are not 
always less energetic than the Z string with the same /3, ^w, Specifically, for n = 1 they 
indeed have lower energy than the corresponding Z strings, but their parameters /5, 6^ lie 
in this case in the non-physical region (see FigJl]), while if they were the true 'W-dressed 
Z strings' they would exist for any parameter values. The special chiral solutions can also 
be found in the physical region (see Fig JTQl) . but only for higher values of n,v starting 
from = 4, z/ = 7 (see FigJl]). In this case they turn out to be more energetic than the 
corresponding Z strings. 

Although the special chiral solutions are not 'dressed' Z strings, one can relate them to 
Z strings very much in the same way as one interpolates between the 'bare' and 'dressed' 
vortices in the Witten model. Given a special chiral solution one can iteratively decrease 
q relaxing at the same time the = condition. This gives solutions with > which 
in the g — limit reduce to the Z string. The current X(g) then shows (see Figs J20!) the 
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FIG. 20: The current X and twist a against q (left) and the magnetic energy versus current (right) 
for the family of the n = v = \ electroweak vortices for (3 = 2 and sin^ 0w = 0.56 that interpolates 
between the Z string for q = and the special chiral solution for q = qi, = 0.15. 

quenched behavior very similar for that for the Witten strings (see Fig l24l) . but it should 
be stressed again that this type of behavior is not generic in the electroweak theory, since it 
requires the fine-tuning of 6^. As is seen in FigEO], the energy of the special chiral solution 
can be higher than for the corresponding Z string. 

D. Infinite Higgs mass limit, /3 — > oo 

In this limit the following constraint is enforced, 

$t$ = /2 + /| = l, (11.43) 

so that one can parametrize the Higgs field amplitudes as 

/i = cose(p), /2 = sine(p). (11.44) 

The equations for the gauge field amplitudes m, Wi, M3, v, f i, are then given by 
fli:T2|) . fliA3D . fHT6D - fHT9D with /i, /a defined by ffTTlill . while the equation for 6 reads 

1 O"^ 1 

- (pQ')' = — (ui cos 26 - M3 sin 2Q)u+ — ivi cos 26 - fa sin 26) v. (11.45) 
p 2 2p^ 
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The constraint (14.201) becomes 



(M1M3 - MgM'J + — {vxv'^ - v^v'^ - g^<$)' = 0. 



11.46) 




1.5 2 

ln(l + p) 



in(l + p) 



FIG. 21: The n = i^ = l, /5 = oo sohitions for cr = 5 (left) and for a = 0.1 (right). The electric 
amplitudes are not shown. 



The boundary conditions for the function G for z/ = n read 



TT 



a5p" + 



2 VP 



+ 



;ii.47) 



which imphes that one always has q = /2(0) = 1. One cannot therefore use q to label the 
solutions, but instead one can use a. The boundary conditions for the gauge field amplitudes 
are given by Eqs. (l5.10l) . (l5.21l) . provided that u = n. For v ^ n the energy density diverges 
at the origin due to the term {y — f3)^/|/p^ in (I4.25p . Therefore, solutions with v ^ n 
become singular for /3 ^ 00, which suggests that they can be viewed as excitations over the 
fundamental regular solution with v = n. 

The typical solutions are shown in Figj2Tl As usual, the limit where a is small corresponds 
to large currents, however, since = 1, the solutions do not show in this limit a region of 
vanishing Higgs field. This can be seen already from Eg. (110.241) . since the size of the $ = 
region scales as ~ X/ a/^ and shrinks to zero when /3 — 00. 

It seems that for (5 = 00 there is no upper bound for cr. Large ex's correspond to small 
currents, the solutions then approaching Z strings. The vacuum angle 7 then tends to zero, 
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while the mass term fl5.22p approaches infinity (since aQ ~ X — and m^j ^ a) in which 
case B changes very fast at small p to approach its asymptotic value. For a = oo one should 
set 0(p) = 0, and then the solutions reduce to Z strings 'in the London limit'. 

XII. SUMMARY AND CONCLUDING REMARKS 

In summary, we have presented new solutions in the bosonic sector of the Weinberg-Salam 
theory which describe straight vortices (strings) carrying a constant electric current. Such 
superconducting vortices contain a regular central core filled with the condensate of massive 
W bosons creating the electric current. The current produces the long-range electromagnetic 
field. The solutions exist for any value of the Higgs boson mass and for (almost) any weak 
mixing angle, in particular for the physical values (3 G (1.5,3.5) and sin^6'w = 0.23. They 
comprise a family that can be labeled by the four parameters in the ansatz (15.61) : o"q, n, u. 
These parameters determine the vortex electric charge density Jq ~ ctq, the electric current 
-^3 ~ CTS; the electromagnetic and Z fiuxes and ^ Zi as well as the vortex momentum 
P ~ cTofTs and angular momentum M ~ (Tq. 

The spacetime vector = (Jq, /s) ~ da is generically spacelike, so that its temporal 
component can be boosted away to give = and o"q = 5^ a. The restframe current I 
can assume any value. In the X — limit the solutions reduce to Z strings, while the twist 
a reduces to the eigenvalue of the linear fiuctuation operator on the Z string background. 

For large currents the solutions show a symmetric phase region of size ~ X where the 
magnetic field is so strong that it drives the Higgs field to zero. However, this does not 
destroy the vector boson superconductivity, since the scalar Higgs field is not the relevant 
order parameter in this case. The current- carrying W-condensate is contained in the very 
centre of the symmetric phase, in a core of size ~ 1/X, while the rest of this region is 
dominated by the massless electromagnetic and Z fields. The symmetric phase is surrounded 
by the 'crust' layer where the Higgs field relaxes to the broken phase while the Z field becomes 
massive and dies away. Outside the crust there remains only the long-range Biot-Savart 
magnetic field, which produces a mild, logarithmic energy divergence at large distances for 
their core. However, finite vortex pieces, as for example vortex loops, will have finite energy. 

Straight, infinite vortices can have finite energy in special cases: either for = 0,7r/2 
when the massless fields decouple, or for ctq = 0, when the current vanishes. The latter 
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case includes Z strings and also the currentless limit of the chiral solutions with (13 = ±o"o- 
However, the chiral solutions are not generic and exist only for special values of 6^. 

The W boson condensate producing the vortex current can be visualized as a superpo- 
sition of two oppositely charged fluids made of and W~ , respectively, flowing in the 
opposite directions. In the vortex restframe the densities of both fluids are the same, which 
is why the total momentum through the vortex cross section vanishes while the current does 
not. Passing to a different Lorentz frame the fluid densities will no longer be equal, which 
will produce a net momentum. The vortex angular momentum can be explained in a similar 
way if one assumes that the fluids perform the spiral motions in the opposite directions. 

Stability of the new solutions is a very important issue. At first one may think that 
they should be unstable, since in the zero current limit they reduce to unstable Z strings. 
However, it seems that current can stabilize them. So far this has been checked only in the 

= 7r/2 limit, where the complete stability analysis has been carried out 22|. It turns out 
then that all negative modes of current-carrying vortices can be removed by imposing the 
periodic boundary conditions, while one cannot do the same in the zero current limit. It 
seems that the same conclusion applies also for 9^ 7^ 7r/2, although the detailed verification 
of this is still in progress. However, the main argument is simple enough and goes as follows. 

As was mentioned above, the Z string perturbation eigenmodes are proportional to 
g±i{(7ot+(T3z) j^ggg Eq.f l8.2p ) with (To = a/o"! — (T^ where cr^ > is the eigenvalue of the spectral 
problem (18. 3p . It follows that all modes with < a are unstable, since ctq is imaginary. 
Since 

27r , 27r , , 

A = — > A^in = — , (12.48) 
0-3 a 

it follows that all unstable modes are longer than Amin- This suggests that one could eliminate 
them by imposing 2;-periodic boundary conditions with a period L < Amin, which would 
leave 'no room' for such modes to exist. However, this would not remove one particular 
homogeneous mode with (T3 = 0, since it is independent of z and so can be considered as 
periodic with any period. 

Now, when passing to current-carrying solutions the situation will be essentially the same. 
It is clear on continuity grounds that, at least for small currents, the solutions will still have 
negative modes with the wavelength bounded from below by a non-zero value Amin- However, 
the crucial point is that the homogeneous mode will disappear from the spectrum - simply 
because the background solutions now have a non-trivial 2;-dependence, implying that all 
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perturbation modes will depend on z as well. This was checked in 22| for = 7r/2 and it 
is very plausible that for generic 6^ the situation will be the same. 

Since the homogeneous mode is absent, imposing periodic boundary conditions with the 
period L = 2TT/a will remove all negative modes. This is somewhat similar to the hy- 
drodynamical Plateau-Rayleigh instability of a water column 39[ or to the gravitational 
Gregory-Laflamme instability of black strings in the theory of gravity in higher dimensions 



40l |. which manifest themselves only starting from a certain minimal length. Since a de- 
creases with current, L then grows. For large currents one has a ~ X^^ (see the caption to 
FigIT]) so that the length of the stable vortex segment scales as 2^ while its thickness grows 
as I (see Eq.f mTm ). 

There could be different ways of imposing periodic boundary conditions on the vortex 
segment. One possibility is to bend it and identify the extremities to make a loop. If the 
vortex carries a momentum P then the loop will have an angular momentum M which 
may balance it against contraction (see Fig l22l) . One may therefore conjecture that such 
electroweak analogues of the 'cosmic vortons' [41] could exist and could perhaps be stable - 
if they are made of stable vortex segments. Of course, verification of this conjecture requires 
serious efforts, since so far vortons have been constructed explicitly only in the global limit 
of the Witten model 42|. However, even a remote possibility to have stable solitons in the 



standard model could be very important, since if electroweak vortons exist, they could be a 
dark matter candidate. 



M 




FIG. 22: Vortex loops balanced by the centrifugal force (left) and vortex segments in a polarized 
medium (right). 



Another possibility to impose periodic boundary conditions is to attach the vortex ends 



68 



to something. It is known that Z strings, since they are not topological, do not have to 
be always infinitely long but can exist in the form of finite segments whose extremities 
look like a monopole-antimonopole pair jl^, This suggests that their current-carrying 
generalizations could perhaps also exist in the form of finite segments joining oppositely 
polarized regions of space, similar to thunderbolts between clouds. If some processes during 
the electroweak phase transition polarize the medium, then the 'vortex thunderbolts' could 
provide a very efficient discharge mechanism - in view of the very large currents they can 
carry. They could be stable in this case and exist as long as the 'clouds' they join are not 
completely discharged. 

There is however an important difference with the ordinary atmospheric thunderbolts, 
which exist only in a polarized medium and emerge when the large electric field between 
the clouds creates the plasma channel and works against the resistance to drive the charges 
through. In superconducting vortices on the other hand the current fiows without any 
resistance and no electric field is needed. In fact, inside the (restframe) vortex the electric 
field is zero. The vortices do not therefore need any polarization of the surrounding medium 
to exist. 





H z y 



H z y 



FIG. 23: Disintegration of an unstable vortex loop (left) and of an isolated vortex segment (right). 



Finite segments of superconducting vortices can probably be created at high temperatures 
or in high energy particle collisions. Once created, the current will start to leave the segment 
through its extremities, so that the latter will emit positive and negative charges (see Fig J2^ . 
Since the current is carried by W bosons, it follows that one vortex end will emit the 
while the other the W~, until the segment radiates away all its energy. In addition, if the 
initial current is large enough, then the vortex will contain a zero Higgs field region, whose 



69 



shrinking will be accompanied by a Higgs boson emission. The vortex segment will therefore 
end up in a blast of radiation, creating two jets of and W~ and a shower of Higgs bosons. 
One can similarly argue that if the vortex segment is created in the form of an unstable loop 
(for example without angular momentum) then it will shrink emitting a shower of neutral 
bosons - Higgs, Z and photons. 

The creation of superconducting vortex segments or loops with their subsequent disin- 
tegration can presumably be observed in the LHC experiments. It is also possible that 
processes of this type could be accompanied by a fermion number non-conservation, since 
already ^ in the zero current limit the vortices (Z strings) admit the sphaleron interpreta- 
tion [13]. However, a special analysis is needed to work out the details of these processes 
and their experimental signatures. We leave this and other possible applications of the 
superconducting vortices for a separate study. 
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Appendix A. Solutions near the symmetry axis 

The boundary conditions at the origin expressed by Eg. (15.91) imply that at small p one 
has 

u = u{0) + 6u, U3 = 1 + 6U3, Ui = 6ui, 
V = 2n — u + 6v, V3 = u + 6v3, vi = 6vi, 

h = 5/1, /2 = /2(0)+5/2, (A.l) 



where the deviations Su, . . . , Sf2 vanish for p ^ and u{0) is a free parameter, while /2(0) is 
also a free parameter if u = n, whereas /2(0) = if z/ 7^ ra. Inserting this to Eqs. (l4.12p - (l4.20p 
and linearising with respect to the deviations the result is as follows. 
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Let us first consider the case where v ^ n and /2(0) = 0. Then the equations read 

ibu)" + - ibu)' = 0, (A.2) 
P 

{6us)" + -^{6u,y = 0, (A.3) 

{6v)" - - {6vy = 0, (A.4) 
P 

{Svs)" - - {6vsy = 0, (A.5) 
P 

1 -n^ 1 

{SfiY + - - - 5/i = -r{cr'{u{0) + ir - (3)6fu (A.6) 
p p'^ 4 

(5/2)" + - (5/2)' - ^—Y^ 5/2 = T(a'(n(0) - 1)^ - /5)5/2, (A.7) 

plus a system of three coupled equations 

{Smf + - {5my = 4 - 5v,), (A.8) 
p p2 

- i (M)' = ^^^(^^^i - i^Su^), (A.9) 

(T2 5u; + -^5t4 = 0. (A. 10) 

P 

From Eqs. (1A.2[) - (]A.5P one finds 

Su = Ai + A2 In p, = 743 + ^4 In p, 

5u = A5 + A6P^ ^^;3 = A + ^8P^ (A.ll) 

while solution of Eqs. flA.6p . flA.7p can be expressed in terms of the Bessel functions, such 
that for small p one has 

= A9p"(l + ...)+Aiop-"(l + ...), 

5/2 = Anp^'-'il + ...) + + •••)• (A.12) 



The three coupled equations (1A.8|) . (]A.9|) . (]A.10P can be reduced to one third order equation 



whose solution can be obtained in terms of the hypergeometric functions. One then obtains 
for small p 



5ui = Ai3P^(l + ...)+^14P"''(l + ...)+^ 



15 



2 2 

5v, = Ai3 p^+'(l + ...) + 77^ ^14 p-"(l + ...) + ^A,,. (A.13) 

1/ + 2 2 — z/ 
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Eqs. ( ]A. 1T]) - (1A. 13p give the most general solution of the equations, since it contains the 



maximal number of the integration constants: Ai, . . . ,Ai^. Now, we have to suppress all 
the solutions that do not vanish for p ^ 0, so that we set 

Ai = A2 = A3 = Ai = A5 = A7 = Aw = A12 = Au = A15 = 0, (A.14) 

assuming that n > u > 0, whereas Aq, Aq, Ag, An, A13 remain free. This gives the local 
solutions in the linear approximation 



u 


= ai, 


Ml 


= as + . . . , 


Us 


= 1, 


V 


= 2n — u + as p^, 


Vl 


_ +2 , 


V3 


= z/ + 04 p^. 


h 


= a5P'' + ---, 




= gp"~" + ... 


containing 6 free parameters: ai = 


m(0), a2 = v4i3, as = Aq, 04 


\i u = n, /2(0) 7^ then the linearized equations for 5u, . 


They split into three coupled groups, whose solutions can be 


gives essentially the same result as 


in (jA.15), 


u = 


ai. 


Ui = 


aa p" + . . . 


U3 = 


1, 


V = 


n + as p^ + . . . , 


Vl = 


p + . . . , 

n + 2 


V3 = 


n + a4 p^ + . . . , 


fl = 


as p" + . . . , 


f2 = 





(A.15) 



(A.16) 



which also contains 6 free parameters: ai, . . . as and q. Eqs. ( 1 A. 15]) and (1A.16P together imply 
Eq. (l5.10p in the main text. 
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These local solutions are obtained in the linearized approximation. However, they can 
be used to take into account all the non-linear terms in the equations when starting the 
numerical integration at p = 0. This is achieved by rewriting the second order field equations 
Eqs. (14. 12l) -f l4. 19p in the first order form, 

y:(p)=^,(p,y,), r,s = l,...16 (A.17) 

where 

y5 = Us, y6 = pu'3, y7 = v, y8 = v'/p, 

1/9 = ^1, yw = v[, yil=V3, y^2 = vyp, (A.18) 

yn = /l/p", 1/14 = /(/P"-' , 1/15 = 1/16 = /2/p'""'''"' 

with the boundary conditions 



1/1 (0) 


= ai, 


1/2(0) = 


0, 


1/3(0) = 


02, 


1/4(0) = z/a2. 


1/5(0) 


= 1, 


1/6(0) = 


0, 


1/7(0) = 


2n — u, 


1/8(0) = 2a3, 


1/9(0) 


= 0, 


l/io(0) = 


= 0, 


l/ii(0) = 


-- 


1/12(0) = 2a4 


1/13(0) 


= ^5, 


1/14(0) = 




1/15(0) = 




1/16(0) = |n - iy\q. 



(A.19) 



The right hand sides J-'r{p,ys) in Eqs.f lA.17l) are defined only for p > 0. For example, one 
has y[ = J-'i = y2/p which is undefined at p = 0. However, using the field equations one can 
check that y[.{0) = for all r. One can therefore set by definition J^r{p = 0, i/s) = 0, and then 
one can integrate Eqs. (lA. 171) with the boundary conditions (1A.19I) starting exactly at p = 0, 
thus avoiding any approximations. Such a procedure is sometimes called desingularization 
of the equations at a singular point. 

Appendix B. Solutions in the asymptotic region 

Solutions of Eqs . (14 . 1 2 1) - (14 . 1 91 ) ■ (14 . 2 1 p have to approach the Biot-Savart configuration 
(I5.17P far away from the vortex core. Therefore one has 

U = U + 6u, U3 = —u + 5%, Ml = 5^1, 
V = C2 + 6v, V3 = -C2 + 6V3, Vi = 6vi, 
fl = 1 + 5/1, /2 = 5/2, (B.l) 
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where 

u = Q\np + ci, (B.2) 

the deviations 6u, ... ,6/2 tend to zero as r — > 00 and ci,C2,Q are the same integration 
constants as in Eq.f l5.17l) . Inserting ( IB.ip to Eqs. (14. 12i) - fl4.20l) and hnearising with respect 
to the deviations, the resulting hnear system reduces to five independent equations plus a 
subsystem of three coupled equations. The first five equations are the following. The linear 
combinations 

6uz = 6u + 6u3, 6vz = 6v + 6V3 (B.3) 

fulfill two independent equations 

{6uz)" + -iSuzy = Isuz, 
P ^ 

{bvz)" --{^vz)' = \5vz. (B.4) 
P ^ 

Their solutions vanishing for r ^ 00 are 

buz = ^ e-™^" + . . . , 6vz = c,^e-^-P + ..., (B.5) 

Vp 

where C3, C4 are integration constants and the dots denote sub leading terms. This corre- 
sponds to the temporal and azimuthal components of the Z field with the mass mz = 1/ V^. 
The linear combinations 

Sua = g^Su - g'^Sus, Sva = g'^Sv - g'^Sv^^ (B.6) 

correspond to the electromagnetic field and satisfy 

{5ua)" + - {Sua)' = 0, 
P 

{5va)" - - {5va)' = 0. (B.7) 
P 

The only solution of these equations that vanishes at infinity is the trivial one, 

Sua = 5va = 0, (B.8) 

since the electromagnetic degrees of freedom are already incorporated into the background 
configuration (15.171) . One has also 

+ - (5/1)' = ^ 5/1, (B.9) 
P ^ 
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whose solution for large p is 



5fi = —e-'^^'P + (B.IO) 



which corresponds to the Higgs boson with the mass = (3/2. 

Let us now consider the equations for 5ui,5vi,5f2 obtained by linearising 
Eqs. f l4. . fl4. 16p . fl4. 181) . The analysis turns out to be more involved in this case. It is 
convenient to replace one of these three equations by the first integral (I4.23p . which gives a 
completely equivalent system. Linearising and defining 

X = y + ^V = m2 (B.U) 



also 



5/2 = (T\X + -y, 5ui = g'^X -uY (B.12) 
2 



the resulting equations read 

X" + (1 + ^) X' - (4 + X) ^ = — y\ (B.13a) 
P X PX 

{Sv,)" - ^ {6v,y -xSvi = cxY, (B.13b) 

^(W-xr'-^X^^, (B.13C) 
P P P 

where C is the same integration constant as in Eq.f l4.23i) . These three equations are equiv- 
alent to one fourth order equation 

o/.2v'^/,^4 6 , 4(c2 - 1)\ ^C2(5 + 4pV^ 



provided that 

S'^.-IP, ^'-X^'-^. (B.15) 

Let us first set C = 0. Then the asymptotic form of the two independent solutions that 
vanish at infinity is given by 

Z= (^Ai/i-5/2(l + ...) + A2Pv^(l + ^ + ...)^exp(- j fidp). (B.16) 

The original field amplitudes can be reconstructed with Eqs. flB.13bl) . flB.13cl) . flB.15p . flB.12p . 

5f2 = -^e-^^''p + ..., 5ui = ^e-^ + Sv^ = c^^pe- ^''p + . . . , (B.17) 

^fp ^/P 
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where Cg, 07,03 are expressed in terms of the two independent integration constants Ai,A2, 
for example Cg = —A2/y/JI. 

Let us now consider the remaining four solutions of Eqs. (1B.13l) - there should be al- 
together six independent solutions, since the equations contain five derivatives and one 
integration constant, C. Two additional solutions are obtained by choosing the plus sign in 
the exponent Eq. (lB.16|) . but they are unbounded and should be excluded from the analysis. 
Yet another solution, always for C = 0, is given by X = 0, 6vi = —C2Y = const., which 
does not satisfy the asymptotic conditions and should be excluded as well. 

The last solution is obtained by setting C ^ 0, which gives rise to a solution that behaves 
a.s Z = Cc2/(/i^Y^) large p, so that 6vi ~ (7 In In p. Let us call this solution C-mode. 
One can think at first that it should also be excluded. However, the constraint equation 
f l4.20p requires that C = 0, and so, even if the C-mode is included, it will be eventually 
removed by the constraint. The constraint is explicitly imposed on local solutions at small 
p, which is sufficient to insure that it holds everywhere, so that it would be redundant to 
impose it again. One can therefore add the C-mode to (IB. 171) . since at the end, when the 
matching is performed to construct the global solutions, it will be removed anyway. 

Since this mode will finally disappear, one may wonder why one should consider it at all. 
However, including it is essential, since it increases the dimension of space of local solutions. 
From the technical viewpoint this mode is not very convenient, though, because it is not 
localized. However, since its main effect is to increase the dimension of function space, any 
other mode with a non-zero overlap with it can actually be employed. Therefore, we simply 



use (IB. 171) with cq, cj, cs considered as three independent parameters, which leads to a good 
numerical convergence. In fact, these coefficients contain p ~ In p, but since this is a slowly 
varying function, one can treat Cq, Cj, Cg as if they were constants. In addition, without 
changing the leading asymptotic terms, one can replace / fi dp in Eq. (]B.17p by f m„ dp. 
Summarizing, the large p behaviour of the solutions is given by 

u = Q\np + ci + ^e-"'^P + ..., = -Q Inp - Ci + ^ e""^^" + . . . , 
V = C2 + c.g'^^pe-'^-P + ..., ^3 = -C2 + c^g^^pe"'^-'' + ..., 

fi = l + ^e-'"H''+..., f, = ^e-^"'-''^ + ..., 

VP VP 

ui = ^e~^"'-''P + ..., Vi = CsVpe-^"'-''' + ..., (B.18) 

VP 
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where Q, Ci, . . . , Cg are 9 independent parameters and m^^ is defined in (IB. lip . This gives 
rise to Eg. (15.2 ip in the main text. 



Appendix C. Superconducting strings in the Witten model 

Numerical construction of Witten's superconducting strings was performed by many au- 
thors [sj] (see for a review). In this Appendix we reproduce these solutions within the 
same approach and using the same notation as in the main text. We tried to make this 
Appendix self-consistent in order to illustrate our procedure on a relatively simple example. 

The U(l)xU(l) Witten model is defined by the Lagrangian |l| 

^^ = -\Y1 Flt^F'-^^^'' + J2 (DM^D^a - U, (C.l) 

a=l,2 a=l,2 

with F^l' = d^A'i^ - dyA^f!:^ and -D^^a = (<9^ - igaAl^^)(j)a where ga are two gauge coupling 
constants. The scalar field potential 

vanishes for |0i| = 1 and 02 = 0, and this minimum is global if > A1A2. The perturbative 
mass spectrum of the field excitations around this vacuum contains two scalar bosons with 
masses ml = Xirjf and m"^ = JTjf— well as a vector boson with the mass = 2g\rH 

and a massless vector boson. Pairs {A^^\(t)i) and (^Jf \ 02) are sometimes called vortex fields 
and condensate fields, respectively. 

Passing to cylindrical coordinates, where x+iy = pe*'^, and making the stationary, axially 
symmetric ansatz 

A(;)dx>' = -{n- v{p)) d^, ct>, = /l(p)e^'^^ (C.3a) 

Afdx^ = - (aodt + a^dz) (1 - u{p)), 02 = /2(p)e^'^"*+^'^^^ (C.3b) 

92 

with n G Z, the field equations 

d^F;^:^ = 2gMz(f>:Dy<f>a), (C.4) 
dU 

D.D'^a = (C.5) 
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reduce to 



^ (p/0' =[y,+ y(/i - ^i) + lf2 ] fu (C.6a) 



^ (p/2)' = (^^^^^ + y(/2^ - Vl) + if,) f2, (C.6b) 

p(^^J = 2glf,v, (C.6c) 

i (/m')' = 2gjflu, (C.6d) 

where = cr'^ — aQ. Depending on sign of cx^, solutions of these equations are called magnetic 
(a^ > 0), electric (cr^ < 0), or chiral (o-^ = 0) 

Equations ( ]C.6P have singular points for p = 0, 00. Constructing their local series solu- 
tions in the vicinity of these points gives the boundary conditions for ^ p —>■ 00, 

1 H — > Ci + Q In p + . . . , 

n + asp^ + ■■■ csv^e"™"'' + . . . , 

aiP+---^/i^l + ^e-'"^^ + ..., 

g + ---^/2^^e-/"^-'^'' + ..., (C.7) 
VP 

where ai, 02, g, Q, ci, . . . C4 are integration constants and 

ml = ml + a\ci + Q\npY. (C.8) 

These boundary conditions imply that the fields are regular at the symemtry axis, while 
at infinity the scalars approach vacuum, the massive vector vanishes, whereas the massless 
vector becomes the Biot-Savart field produced by the string current. The current density is 

J. = = 2^23?(#;/^.02), (C.9) 

integrating which over the x, y plane gives 

[ J^d^x = -'^^^, (C.IO) 

J 92 

with a = 0, 3. Here Jq is the charge per unit vortex length while Is is the total current 
through the vortex cross section. As Eq. llC.SP shows, the long-range Biot-Savart field affects 
the mass of the second scalar. For o"^ > this only improves the field localization, but 
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for 0"^ < the term becomes negative at large distances, making the condensate field 
oscillate at infinity, which renders the energy strongly divergent. Such solutions are more 
difficult to handle, and so we shall concentrate on the magnetic and chiral cases. Since 



the structure of the fields ( IC.SP is invariant with respect to Lorentz boosts along the z- 
axis, this can be used in the magnetic case to pass to the rest frame where la = with 
T = —2T[Qa/g2 where the 'twist' a = \fo^. 

To solve equations flC.6|) the relaxation method is usually employed j^, Q]. Within this 
method it is very natural to find first of all the 'dressed' currentless solution with -u = cr^ = 0, 
/2 7^ 0. After this one can increase to obtain current-carrying strings with m 7^ 0. 

However, we use instead the multiple shooting method in which the local solutions in 
( ]C.7[) are numerically extended to an intermediate region where the matching conditions are 
imposed. To fulfil the 8 matching conditions for the 4 functions and their first derivatives 
we have on our disposal 9 free parameters: 8 integration constants in (1C.7|) and also a. This 
leaves after the matching one free parameter to label the global solutions in the interval 
p G [0, 00). We choose this parameter to be the value of the condensate scalar at the origin, 
(I = /2(0). Within our approach it is then natural to start not from the 'dressed' solution 
that is apriori unknown, but from the known ANO vortex. 

We therefore set g = and then one has /2 = 0, m = 1 everywhere so that the condensate 
fields vanish. The remaining equations reduce to the ANO system whose solution is the 
ANO vortex 21]. For g ^ 1 one can expect the solution to be its small deformation that 
can be considered within the linear perturbation theory. Linearising the equations with 
respect to the deformations the ANO background, one can see that it is consistent to set 
the deformations of u,v, fi to zero, while the deformation of /2 satisfies 



^ W2) = (^^ -f^2+7//j5/2. (C.ll) 

This can be viewed as the spectral problem with the eigenvalue o"^. This problems admits 
a bound state with o"^ > 0, which leads to the dispersion relation 



This interpretation of this is two- fold. If the wave- number is small, (j| < cr^, then the 
requency (Tq is imaginary, implying that the ANO vortex within Witten's model is unstable 
l|. On the other hand, for (t| > cr^ the frequency is real and the solution of (1C.11|) describes 

a small stationary deformation of the ANO vortex by the condensate. 
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The perturbative bound state can be used as the starting configuration to iteratively 
construct the solution of the full non-linear system, still for small q. After this one can 
iteratively increase g, which gives the generic solutions. The typical solution is shown in 
FigEH One finds then that the restframe current X first increases with q but then starts to 
quench, and finally both X and a (but not Q) vanish for some as shown in FiglJH 




0.5 1 1.5 2 2.5 3 0.05 0.1 0.15 0.2 0.25 0.3 



ln(l + p) q 

FIG. 24: Left: profiles of a typical current-carrying solution of Eqs. ()C.6P for Ai = 0.1, A2 = 10, 
r/i = 1, = 0.31, 51 = (72 = 1) 7 = 0.6 and a = 0.12. Right: the restframe current X and the twist 
a against the parameter q for the same parameter values as in the left panel. 

It is worth reminding that in the Weinberg-Salam theory we had a completely different 
picture, as shown in Figj6]- when q approaches its maximal value the current 1{q) does 
not quench but bends back to smaller values of q and continues to grow. 

Getting back to Witten's model, the solution for q = q-k (it looks qualitatively similar to 
other solutions) has = and so it is chiral, with values of ctq = ±(73 related to each other 
by Lorentz boosts. In this case one has X = but the current does not actually vanish 
and becomes an isotropic vector. The particular member of the chiral family is obtained for 
c"o = c^s = 0, in which case one has Jq, = so that the solution is truly currentless, although 
it contains a non-trivial condensate. This is the 'dressed' string. 

Summarizing, when q increases, the Witten strings interpolate between the 'bare' ANO 
vortex for g = and the 'dressed'/chiral vortex for g = g^. Since the restframe current 
vanishes at the ends of the interval [0,g^], it passes through a maximum somewhere in 
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between and so there is an upper bound for it. The magnetic energy, which coincides with 
the total energy for the currentless strings, decreases along the interpolating sequence, as 
shown in Figl25l The 'dressed' vortex is therefore less energetic than the 'bare' one. 





FIG. 25: Magnetic energy (left) and the restframe current X and q = /2(0) (right) against the twist 
a for the same parameter values as in Fig[ 



Within the standard approach based on the relaxation method the same solutions are 
obtained in the opposite order. First one finds the 'dressed' solution with g = g^, a = 0, 
and then one increases the twist a. The functions g(cr), X((t) then exhibit the characteristic 
behaviour shown in Figl25] (which agrees with Fig. 5. 3 in Ref.[5j] where a is called 'winding 
number density'). The condensate then always decreases, while the current passes through 
a maximum whose value can be related to the parameters of the 'dressed' solution 5|. 

This scheme does not apply in the Weinberg-Salam theory, where the 'dressed' currentless 
solutions do not generically exist. We are therefore bound to start from the embedded ANO 
vortex with g = 0. Increasing q we find that its maximum value exists, although it 
corresponds not to the 'dressed' currentless solution but to a turning point after which g 
starts decreasing again, while the current continues to grow. 
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